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There  has  been  much  interest  in  the  study  of  the 

properties  of  electrolytes  in  dipolar  aprotic  solvents.   Two 

such  solvents  were  chosen  for  this  study:   sulfolane,  and 

l-methyl-2-pyrrolidinone  (NMgPy) .   The  electrolyte  chosen  for 

this  study  was  hydrogen  chloride  (HC1),  a  stong  electrolyte 

in  water.   The  conductimetric  behavior  of  HC1  was  studied 

in  sulfolane  and  sulfolane-water  solvents  at  25,  30,  and 

40°C,  and  in  NMgPy  and  NM2Py-water  solvents  at  25°C.   A 

study  was  made  of  the  conductimetric  behavior  of  sodium 

chloride  (NaCl),  sodium  acetate  (NaAc) ,  tris(hydroxy- 

methyl  aminomethane  hydrochloride  (Tris.HCl),  acetic  acid 

(HAc) ,  and  HC1  in  the  mixed  solvent  containing  80  weight-% 

2-methoxyethanol  and  20  weight-%  water  (80  W%  MetOH-HgO). 
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Conductances  were  measured  with  fill  cells  of  the  Kraus 
type,  in  conjunction  with  a  transformer  bridge  based  on  the 
design  of  Janz  and  Mclntyre.   Measurements  were  made  at 
2.5,  5.0,  and  10.0  kHz  and  extrapolated  to  obtain  the  true 
ohmic  resistance  by  the  method  of  Hoover.   The  conductance 
data  were  analyzed  with  the  full  Pitts  equation,  setting  the 
ion  size  parameter  equal  to  the  Bjerrum  critical  distance. 
Analyses  were  also  made  with  log-log  plots  and  the 
Shedlovsky-Kay  extrapolation.   A  comparison  of  the  full  and 
expanded  conductance  equations  was  also  made. 

The  study  of  the  sulfolane  solvent  systems  showed 
increasing  ion  association  and  decreasing  limiting  equi- 
valent conductance  with  increasing  sulfolane  content .   In 
pure  sulfolane,  HC1  was  found  to  exist  primarily  in  an 
associated  form,  with  evidence  of  homoconjugation  with 
chloride.   There  was  evidence  of  a  change  in  the  dis- 
sociation process  for  HC1  in  pure  sulfolane  with  increasing 
temperature . 

The  study  of  HC1  in  the  NMgPy  solvent  systems  showed 
a  sharp  decrease  in  limiting  equivalent  conductance  with 
increasing  NM„Py  content ,  coupled  with  a  gradual  increase 
in  ion  association.   As  opposed  to  sulfolane,  NM„Py  gave 
evidence  of  being  able  to  solvate  the  proton.   For  both 
sulfolane  and  NM„Py,  the  electrostatic  models  of  ion 
association  were  generally  unable  to  predict  the  amount  of 
association  that  was  found. 
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Both  of  the  dipolar  aprotic  solvents  studied  gave 
evidence  of  disrupting  the  water  structure,  with  sulfolane 
causing  a  decrease  in  the  over-all  structuring  of  the 
solvent.   The  NM„Py  was  found  to  cause  an  increase  in  the 
structuring  of  the  solvent ,  with  evidence  of  the  formation 
of  an  associated  species  of  the  form  NM^Py^HgO. 

In  the  study  of  the  80  W%  MetOH-H20,  all  electrolytes 
studied,  except  HAc ,  were  found  to  be  slightly  associated. 
The  conductance  of  Tris-HCl  in  water  was  studied,  and  a 
comparison  between  the  two  systems  was  made.   Acetic  acid 
was  found  to  exist  in  a  primarily  associated  form,  with  the 
dissociation  constant  comparable  to  that  found  in  other 
such  alcohol-water  solvents. 
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CHAPTER  I 
INTRODUCTION 

The  last  two  decades  have  seen  a  growing  interest  in 
that  class  of  solvents  characterized  as  being  both  dipolar 

and  aprotic.   Some  of  these  solvents,  such  as  acetonitrile , 

1-6 

have  been  extensively  studied;     others  have  been  less  well 

characterized.   One  such  is  sulfolane  ( tetrahydrothiophene- 
1-1-dioxide,  tetramethylene  sulfone).   One  reason  for  this 
lack  of  research  may  be  that  sulfolane  is  a  solid  at  room 
temperature,  requiring  measurements  at  temperatures  of 

30°C  or  above.   Yet  despite  this  drawback,  sulfolane  remains 

7—8 
an  excellent  solvent.      Possessing  a  moderately  high 

9 
dielectric  constant,   sulfolane  would  be  expected  to  reduce 

the  ion-ion  interactions  found  in  solvents  with  smaller 

dielectric  constants.   With  its  high  dipole  moment    it 

would  be  expected  to  increase  ion-dipole  interactions.   For 

both  these  reasons,  sulfolane  is  expected  to  have  a  moderate 

11-12 
solvating  power.   Sulfolane  has  been  characterized     '  as 

having  extremely  weak  acid-base  properties,  facilitating  the 
study  of  strong  acids  and  bases. 

Another  of  the  lesser  dipolar  aprotic  solvents  is 
l-methyl-2-pyrrolidinone.   For  the  same  reasons  cited  for 
sulfolane  it  is  expected  to  have  moderate  solvating  power; 

-1- 


-2- 

Q 

it  is,  however,  more  basic  than  sulfolane,   and  has  been 

13 

shown  to  enter  into  extensive  hydrogen  bonding  with  water, 

14 
something  which  is  not  true  of  sulfolane.     It  would  be  of 

interest  then  to  compare  the  behavior  of  electrolytes  in 

these  two  dipolar  aprotic  solvents. 

A  far  more  common  class  of  solvents  is  that  considered 

to  be  protic,  the  most  common  of  these  being  water.   One  of 

the  lesser  studied  of  this  class  of  solvents  is  the  mixed 

solvent  consisting  of  80  weight  percent  2-methoxyethanol 

(methyl  cellosolve)  and  20  weight  percent  water.   First 

15 

studied  by  Simon  "    in  the  middle  fifties,  a  number  of 

studies  have  since  been  published  regarding  this  solvent, 

but  questions  about  the  behavior  of  electrolytes  within  it 

17 
remain. 

19 
Conductance  studies  have  been  shown    to  be  an  excellent 

tool  for  characterizing  the  behavior  of  electrolyte  solu- 

20 
tions,  due  in  part  to  the  great  accuracy  attainable.     It 

is  the  major  tool  employed  in  these  studies. 


CHAPTER  II 
THEORY 

General  Background 
The  resistance  of  a  uniform  cross  section  of  a 
homogeneous  substance,  with  length  1  and  area  A  is  given 
by  the  equation: 

R  =  AT  C1) 

where  k  is  a  constant  called  the  specific  conductance 
and  is  characteristic  of  the  material  under  study.   If  a 
conductivity  cell  had  a  uniform  cross  sectional  area  and 
electrode  separation,  then  by  measuring  the  resistance  of  a 
solution  we  could  obtain  the  specific  conductance.   Since 
cell  geometry  varies  with  manufacture,  however,  it  becomes 
necessary  to  introduce  a  parameter  called  the  cell  constant. 
This  will  have  units  of  length  divided  by  area  and  can  be 
determined  by  measuring  the  resistance  of  a  solution  whose 
specific  conductance  is  known.   The  relationship  between 
specific  conductance  and  cell  constant  is: 

<=|  (2) 

Since  the  conduction  of  electricity  by  an  electrolyte 
solution  depends  on  ions  as  the  charge  carriers,  the 
specific  conductance  will  vary  with  the  concentration  of  ions 
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By  dividing  the  specific  conductance  by  the  concentration 
we  obtain  a  conductance  per  mole,  useful  for  making  compar- 
isons between  solutions  of  differing  concentrations.   This 
is  the  molar  conductance,  also  known  as  the  equivalent 
conductance.   It  is  defined  by  the  equation: 

A  =  looo,  =  loooe 

c       Re  v 

where  c  is  the  molar  concentration.   The  resistance  is 
usually  corrected  for  the  conductance  due  to  the  solvent: 

1=1 L_  (4) 

R    R     R   ,  y     } 

m    solv 

where  R   is  the  measured  resistance  and  R   ,   is  the 
m  solv 

solvent  resistance.   In  this  study  the  terms  equivalent 
conductance  and  molar  conductance  are  used  interchangeably, 
since  for  a  1:1  electrolyte,  such  as  is  studied  here, 
equivalents  per  liter  and  moles  per  liter  are  identical. 

It  was  in  1883  that  Arrhenius,  by  studying  the  varying 
conducting  power  of  solutions  came  to  the  conclusion  that 
electrolytes  could  supply  varying  amounts  of  ions  depending 
on  the  electrolyte  concentration.   He  concluded  that  as  the 
concentration  tended  toward  zero  (infinite  dilution),  the 
degree  of  dissociation,  a,  tended  toward  unity.   Noting  that 
at  a  fixed  temperature  A  increased  as  concentration 
decreased,  he  ascribed  this  to  increasing  dissociation.   Thus 
at  infinite  dilution,  A  would  reach  a  maximum  value,  A  , 
which  would  correspond  to  the  conductance  for  a  totally 
dissociated  electrolyte,  and  the  degree  of  dissociation 
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at  any  given  concentration  could  be  expressed  as  a  ratio 
of  the  molar  conductance  at  the  concentration  to  the  molar 
conductance  at  infinite  dilution: 

a  =  —  (5) 

A° 

From  this,  Ostwald  later  derived  the  famous  dilution 

law : 

2  .2 

v      =        a   c    -     A  c 
d  "  CI  -  «)  "  Ao(Ao  _  A)  <6> 

This  was  in  accord  with  the  Arrhenius  proposition  that 

at  low  concentrations  the  molar  conductance  should  vary  in 

21 
a  linear  manner  with  the  concentration." 

However,  Kohlrausch,  a  contemporary  of  Arrhenius, 

had  made  some  precise  measurements  of  the  conductance  of 

solutions  at  low  concentrations,  and  he  came  to  the  conclu- 

i  <->2 

sion  that  A  varied  with  c2  and  followed  the  relation:" 

A  =  A°  -  Ac*  (7) 

There  were  sets  of  data  which  conformed  to  one  theory; 
there  were  data  sets  which  conformed  to  the  other.   This 
seeming  paradox  held  sway  for  many  years. 
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Measurement  of  Conductance 
Before  proceeding  further,  a  few  words  are  in  order 
on  how  the  conductance  of  a  solution  is  measured.   Such 
measurements  are  extremely  sensitive  and  require  very  care- 
ful technique  in  order  to  insure  accurate  results.   As 

23 
Shedlovsky  has  pointed  out,"   the  resistance  measured  for 

a  solution  will  change  approximately  2%   per  °C;  to  have  a 
precision  of  0.01%  thus  requires  that  the  temperature  be 
regulated  to  0.005°C  or  better.   This  is  but  one  of  several 
factors  that  need  to  be  considered  in  conductance  measure- 
ments . 

When  determining  conductance  one  is  actually  measuring 
the  resistance  of  a  solution.   As  with  most  precision 
resistance  measurements,  some  variant  of  Wheatstone's  bridge 
is  usually  used.   A  schematic  diagram  of  such  a  bridge 
appears  in  Figure  1.   The  circuit  is  such  that,  when  the 
detector  registers  a  null  (i.e.,  there  is  no  current  flowing 
in  the  circuit)  the  following  balance  condition  is  met: 

R    R 

R..  and  R„  are  precision  resistors  whose  values  are  known 
(and  preferably  equal);  thus  at  null,  the  cell  resistance 
is  equal  to  the  resistance  set  on  the  standard. 

Such  a  bridge  was  designed  primarily  for  direct  current 
measurements.   However,  if  a  direct  current  is  applied  to 
an  electrolyte  solution  there  will  be  migration  of  ions  to 
the  electrodes;  this  electrolysis  will  produce  a  change  in 


-7- 


Figure  1 
Wheatstone's  Bridge 
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Figure  2 
Jones  —  Joseph  Conductivity  Bridge 
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the  solution  concentration  and  a  concomitant  change  in 
the  resistance. 

Initially,  it  was  Kohlrausch  who  proposed  to  solve  this 
problem  by  the  use  of  alternating  current  for  resistance 
measurements.   This,  however,  introduces  a  problem  of  its  own. 
Given  the  parallel  plate  configuration  of  most  conductance 
cells,  there  will  arise  a  capacitance  due  to  the  flow  of 
alternating  current  through  the  cell.   This  means  that  the 
bridge  must  not  only  balance  resistances,  it  must  now  balance 
capacitances  as  well.   As  a  result,  modern  bridges  are 
designed  to  balance  impedance. 

A  great  deal  of  the  early  work  on  bridge  design  and 

20         24 
conductance  measurement  was  done  by  Shedlovsky,    Parker, 

Jones  and  his  students  °'   '   '   '    and  others.   Much  of 

30 
this  early  work  has  been  summarized  elsewhere    and  will  not 

be  covered  in  any  great  depth  here.   A  typical  bridge  circuit 

95 

designed  to  measure  the  impedance  of  a  conductance  cell    is 

presented  in  Figure  2.   In  more  recent  times,  Ives  and  his 

31 
coworkers  have  demonstrated    the  usefulness  of  the  trans- 
former bridge  circuits  found  in  commercial  impedance 

comparators.   This  study  utilized  just  such  a  bridge, 

32 

designed  by  Janz  and  Mclntyre.     It  will  be  described  in 

more  detail  in  a  later  chapter. 

When  measuring  resistance  with  such  a  bridge,  a 
dependence  of  resistance  on  the  frequency  of  the  alternating 
current  used  is  unavoidable.   For  precise  measurements  to 
have  any  meaning,  there  must  be  some  way  to  eliminate  this 
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frequency  dependence  and  obtain  the  true  ohmic  resistance. 
A  large  amount  of  the  frequency  dependence  has  been  found 
to  be  caused  by  extraneous  factors  and  can  be  removed  by 
proper  cell  design   ''    and  the  use  of  a  proper  thermo- 
stating  medium. 

Platinization  of  the  electrodes  helps  to  reduce  the 

34 
remaining  frequency  dependence  greatly,    and  is  of  great 

aid  when  measuring  the  conductance  of  concentrated  solutions; 

in  dilute  solutions,  the  large  surface  area  of  the  platinum 

black  favors  adsorption  of  the  electrolyte  under  study. 

Consequently,  in  this  study,  shiny  platinum  electrodes  were 

used. 

A  conductance  cell  may  be  schematically  represented  as 

shown  in  Figure  3: 

R„ 

Rl 


1 


C„ 


Figure  3 
Circuit  equivalent  to  conductance  cell 


3-1 
This  circuit  was  proposed  by  Ives  and  his  coworkers'  "  and 

endorsed  by  Robinson  and  Stokes,  who  added  the  Warburg 

34 
impedance,  -W-.     In  this  circuit,  R..  is  the  true  ohmic 

resistance  of  the  cell,  which  is  what  we  wish  to  determine. 

This  is  independent  of  frequency  over  the  limited  range  of 

measurement  (between  1  and  10  kHz).   In  series  with  R.. 
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is  the  capacitance  C.  ,  due  to  the  double  layer.   This  is 
also  expected  to  be  independent  of  frequency  in  the  range 
of  frequencies  used  in  this  study.   It  is  this  capacitance 
which  in  fact  is  responsible  for  transporting  through  the 
double  layer  the  current  passing  through  R- ;  it  has  been 
shown  that  a  cell  has  a  definite  resistance  with  only  a  few 
millivolts  of  potential  applied,  too  small  a  potential  to 
cause  electrolysis  at  the  electrodes.   There  is,  however,  a 
small  amount  of  electrolysis  which  takes  place,  probably  due 
to  discharge  of  solvent  ions,  and  this  is  represented  by 
"faradaic  leakage"  in  parallel  to  the  double  layer.   It  is 
represented  by  resistance  R2  and  the  Warburg  impedance,  both 
of  which  are  expected  to  vary  with  u  2 .   Capacitance  C„ 

is  the  capacitance  due  to  the  leads  and  connections,  and  in 

34 
cells  such  as  were  used  in  this  study  it  is  negligible. 

31 
Ives    has  assumed  a  negligible  Warburg  impedance  and 

determined  the  frequency  dependence  of  the  measured 

resistance  to  be  of  the  form: 

Rm  =   Rl   +  — -2  (9) 

1   +   aoo 

where : 

R1R2   Cl 
Rl   +   R2 

and : 

to    =    2ttv  ,    v    being   the    frequency. 
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35 
Hoover    has  compared  a  number  of  models  for  bridge 

balance  conditions  and  has  come  to  the  conclusion  the  model 

proposed  by  Ives  is  of  sufficient  accuracy  to  use  for  the 

extrapolation  to  remove  the  frequency  dependence  and  obtain 

the  true  ohmic  resistance.   A  more  detailed  discussion  of 

this  frequency  extrapolation  method  appears  in  a  later 

chapter. 

Over  the  years  there  have  been  a  large  number  of  cell 

designs  for  the  measurement  of  electrolytic  conductance.   For 

36 
the  purposes  of  this  study  a  cell  of  the  Kraus  design    was 

used.   This  type  of  cell  helps  to  eliminate  several  sources 

of  error,  including  some  sources  of  frequency  dependence 

in  the  measured  resistance  alluded  to  above.   A  discussion 

of  this  style  of  cell  appears  in  the  experimental  section 

along  with  a  photograph  of  the  cell  actually  used. 
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Conductivity  Equations 
As  mentioned  previously,  in  the  early  days  of 
conductance  studies  there  was  something  of  a  paradox.   The 
Ostwald  dilution  law,  based  on  the  Arrhenius  theory  of  dis- 
sociation, indicated  a  linear  dependence  of  equivalent 
conductance  on  concentration  and  there  were  data  to  back  it 

up.   Yet  Kohlrausch  showed  a  dependence,  not  on  c  but  on 

i 
c2 ,  and  had  data  to  confirm  the  variation.   Kohlrausch  had 

also  proposed  his  law  of  independent  migration,  which, 

succinctly  stated,  says  that  the  limiting  conductance  (the 

conductance  at  infinite  dilution)  is  equal  to  the  sums  of 

the  limiting  conductances  for  the  individual  ions: 

A°  =  A°  +  A°  (10) 

there  being  no  interaction  between  moving  ions  at  infinite 
dilution.   Since  it  is  unlikely  that  both  ions  will  contri- 
bute identically  to  the  total  conductance,  the  fraction 
carried  by  one  ion  may  be  expressed  by  the  transport  number 
t°: 

*?-;!'  (11) 

The  expression  for  the  negative  ion  is  identical,  and  thus 
from  Equation  10  one  finds: 

t°  +  t°  =  1  (12) 

The  paradoxical  problem  continued  to  plague  conduc- 
tivity through  the  first  several  decades  of  this  century. 
In  1923,  Debye  and  Hiickel    proposed  their  theoretical 
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treatment  of  interionic  interactions,  and  this  was  further 

33 
applied  to  conductance  studies  in  1927  by  Onsager.     This 

combined  Debye-Hiickel-Onsager  theory  was  the  first  major 

attempt  to  theoretically  account  for  the  variation  of 

conductance  with  concentration. 

In  an  electrolyte  solution  placed  under  the  influence 
of  an  external  electric  field  there  are  two  major  interionic 
interactions:   the  electrophoretic  effect  and  the  relaxation 
effect.   These  have  been  extensively  discussed  elsewhere   ' 
and  thus  will  be  dealt  with  in  a  more  succinct  manner  here. 

The  electrophoretic  effect  arises  from  viscous  drag. 
An  ion  moving  through  a  solution  will  tend  to  drag  along 
solvent  molecules  with  it,  both  through  frictional  forces  and 
through  having  solvent  molecules  associated  with  the  ion 
(solvating  the  ion).   Neighboring  ions  will  thus  experience 
a  net  solvent  flow,  in  the  direction  of  their  movement  if 
they  are  the  same  charge  type  as  the  central  ion,  and 
against  their  motion  if  of  opposite  sign  to  the  central  ion. 
Quite  clearly,  this  effect  will  decrease  as  the  concen- 
tration of  ions  decreases,  and  will  reach  zero  at  infinite 

dilution. 

34) 
It  has  been  shown    that  the  electrophoretic 

increment  to  the  velocity  of  the  ion  is  given  by  the 

equation : 

n,  n-1,      n-1.  . 
ro       z  (z    k   -  z    k  ) 

v   =  I      A   •  -=— = - —  (13) 

-i   n        n/        \ 
n=l  a    (z_  -  z  ) 
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where  n  arises  from  Boltzmann  distribution  functions,  k 
is  the  force  acting  on  the  ion,  and  A   is  given  by: 

A   =  ^i^-l^L)*'1®    (««0  (14) 

n    n  !6iTri  "-ekT-'     nv  '  v  J 

The  function  $  (Ka)  is  given  by: 


$  (<a)  =  (Ka)zf^- )nS  (Ka) 

nv  J         v  '     y\.   +  Ka;   nv   y 

(15) 

„  „  -nKr 

o  /   \     n-2  re      ^ 

S  (Ka)  =  a      ^ —  dr 

nv  '  I      n-1 

a  r 

k  in  this  case  is  not  the  specific  conductance,  but  the 

Debye-Hiickel  reciprocal  radius  of  the  ionic  atmosphere,  and 

a  is  the  distance  of  closest  approach  of  the  ions;   k  is 

given  by: 

K  =    f8irNe  I  yj  (16) 

K         ^lOOOekT-'  {lb) 

with  the  remainder  of  the  symbols  having  their  usual 

meanings. 

The  relaxation  effect  arises  as  follows.   The  ion 

being  studied  is  at  the  center  of  an  atmosphere  of  spherical 

distribution.   If  the  central  ion  moves,  the  atmosphere 

will  move  with  it,  but  a  little  more  slowly,   As  a  result, 

the  ion  feels  a  net  restoring  force  in  the  direction 

opposite  to  its  motion,  until  the  atmosphere  has  once  again 

assumed  the  spherical  symmetry .   If  the  field  acting  on  the 

ion  is  designated  by  X,  then  the  relaxation  field  may  be 

denoted  as  AX.   Initially,  Debye  and  Hiickel  examined  the 

relaxation  field,  but  it  was  elaborated  by  Onsager.   For 
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the  1:1  electrolytes  we  are  concerned  with  in  this  study, 
the  expression  derived  by  Onsager  is : 

AX  _  z  e       qK  M  . 

X  "'  3ekT  '  1  +  /q  ( -1  ° 

where  q  is  a  function  of  charge  z  and  ionic  conductance 

at  infinite  dilution  \    .      For  a  1:1  electrolyte  its  value 

is  §.   The  velocity  of  the  ion,  corrected  for  the  relaxation 

effect  becomes : 

2 


z.|  z^e 
x   3ekT       ~75- 


"<   =  v°  (i  +  4Sr-  ■   r^xr)  (W) 


The  effective  force  acting  on  a  moving  ion  must  be 

40 
equal  and  opposite  to  the  force  of  the  external  field    and 

may  be  expressed  as : 

-z.eX  =  FV1S  +  Fl  +  F"  (19) 

3  3  3  3 

where  F.    is  the  viscous  drag  of  the  solvent  on  the  moving 

ion,  F.  is  due  to  the  electrophoretic  effect,  and  F.  is  due 
J  3 

to  the  relaxation  effect.   It  might  be  noted  that  at 

infinite  dilution  F.  and  F '.  decrease  to  zero,  and  it  is  thus 

J       3 

F.    which  determines  the  limiting  conductance,  A... 
J  .J 

Combining  all  of  the  above  factors,  Onsager  derived 
what  has  come  to  be  known  as  the  Debye-Huckel-Onsager 
limiting  law: 

2    2  .o  v2 

A    =    A°    -   ^=    •     /^      -   7^f(2z)K  (20) 

3ekT         1    +    /q         OtthN 

For  a  1:1  electrolyte,  this  reduces  to: 

A  =  A°  -  Sc*  (21) 
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where : 


and 


S  =  aA°  +  6 

4     z3 
a  =  82.0460  x  10   =^— =  (22) 

(eT)1'13 


3 
6  =  82.487  — - r 


n(eT)2 
As  can  be  seen,  the  final  equation  has  the  same  form  that 

was  empirically  determined  by  Kohlrausch,  the  same  depen- 

x 
dence  on  c2.   But  now,  the  constant  S  has  a  theoretical 

basis . 

Since  the  equation  was  first  proposed  by  Onsager,  there 

have  been  numerous  treatments  devised.   Most  of  them 

concentrate  on  elaborating  the  terms  for  the  relaxation 

field,  while  the  electrophoretic  effect  has  also  been 

20 
modified.   Shedlovsky""1   proposed  a  form  of  the  limiting  law 

suitable  for  extrapolation  to  determine  A  .   Robinson  and 

41 
Stokes    did  much  the  same  but  included  a  parameter  to 

42 
reflect  finite  ion  size.   Shedlovsky  and  Kay  devised    an 

extrapolation  for  very  weak  (i.e.,  highly  associated) 

electrolytes,  based  on  an  earlier  treatment  by  Fuoss  and 

43 
Shedlovsky . 

Today,  there  are  two  main  equation  systems  for 

analyzing  conductance  data.   One  is  the  result  of  many 

44  ^5_47 
successive  treatments  by  Fuoss  and  his  coworkers. 

48-50 
The  other  is  based  on  a  theoretical  treatment  by  Pitts. 
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Both  theories  are  quite  complex,  with  the  main  equations 

containing  several  exponential  integrals.   Both  have 

51  5° 
been  expanded   '  "  to  the  same  form: 

A  =  A°  -  Sc~  +  Ec  log  c  +  J  c  -  J2c1-5  (23) 

In  both  terms,  the  constants  S  (the  Onsager  limiting  law 
slope)  and  E  are  identical.   The  terms  J-  and  J^  differ 
depending  on  the  theory  used;  it  is  these  terms  which  contain 
the  so-called  ion-size  parameter.   This  will  be  discussed 

in  depth  in  the  discussion  chapter. 

40 
The  equations,  in  the  expanded  format,  are  as  follows: 

2  2  2 

ba  =  |^-  =  16.7099  x  104  ~ 

k   =    50.2916  — ^ 


(eT)2 


E  =  EXA0  -  E2 


E-.  =  (lC»^   =  2.94257  x  1012  5__. 

1      24c  (eT)3 


r-  5 

E   =  KabS  =  4> 33244  x  iq7 


16c2  n(eT)^ 


J   =  2E1A°  (ln(^J)  +  Ax)  +  2E2  (A   -  ln(^x)) 

c2  c2 

j2  =  i££b  (4EiAoA3  +  2E^    _  ^ 


c 

The  delta  functions  are  complex,  and  may  be  obtained  in 
Reference  40.   Their  values  depend  on  the  theory  being  used 
They  are  for  the  Pitts  equation  and  the  latest  of  the  Fuoss 
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44 
equations,  that  of  Fuoss  and  Hsia,    both  of  which  have 

been  expanded  by  Fernandez-Prini  and  Prue. 

In  recent  years,  there  has  been  a  great  deal  of 
discussion  as  to  which  of  the  theories  and  its  associated 
equations  best  fits  the  experimental  data  available.   What 
is  truly  remarkable  is  that  this  has  been  going  on  for 
better  than  a  decade  and  a  half  and  there  is  still  no 
conclusive  proof  that  one  theory  is  better  than  the  other. 
That  being  the  case,  the  theory  used  for  data  analysis  is 
really  serendipitous,  depending  on  individual  whim. 

As  will  be  mentioned  in  a  following  chapter,  the 
expanded  equations  given  above  are  approximations  only  and 
tend  to  break  down  at  high  concentrations  (i.e. ,  where  Ka 
becomes  relatively  large).   For  this  reason,  the  data  of 
this  study  were  analyzed  with  the  full  Pitts  equation 
described  below. 

The  Pitts  conductivity  equation  takes  the  general 
form : 

A/c 

y)(/2 

(24) 


A   =    A°    -    A°    ( ^- +    BcSj    +    G/c"    (M    I      ,    - 

+    y)(/2    +    y)  •      lJ  H1    +    v) 

A/c  H/c 


where 


(1+   y)2(/J+   y)    ~    <1+   *>    Tl) 

A    =    H(/2    -    1) 

B    =    3H2 

ekTHNIO9         „    .       ..  ,        .,  .    , .    . . 
~ —    ,    C   is   the   velocity   of    light 

TTnCT 


-19- 


2  2 

„    z  e  k 
ti.    — 


3ekT/c 


9     9/2  -  10  +  y(3/2  +  1)  +  2y2       eY       , 
Sl_      8(l+y)2(/2+y)2       -4(1  +7)  Ei(y) 


+ 


TH  = 


1    4(/2  +  y 


(1  +  y)  (/2  +  y) 

/2  ,    e(1  +  3)y 
+  16  ((1  +  y)(/2  +  y))  Ei  l<&    +    1>y) 

-  +  |  (3^e<g  +  1)y  E   ((p  +  1)y)_  2eyE.(y)) 


E.(x)  =  / du 

x 

2  2 
y  =  ica,  where  k  =  (-10ooekT  J 

»"* 

All  other  symbols  have  their  usual  meaning.   Values  for  the 
physical  constants  are  taken  from  Reference  53. 

This  equation  applied  to  dissociated  electrolytes.   To 
make  it  applicable  to  associated  electrolytes,  concen- 
trations (c)  were  replaced  by  ionic  concentrations  (etc), 
where  a  is  the  degree  of  dissociation.   Similarly,  the 
calculated  equivalent  conductance  calculated  by  the  program, 
A.,  which  is  the  conductance  for  a  completely  dissociated 
electrolyte,  was  multiplied  by  the  degree  of  dissociation 
to  get  the  equivalent  conductance,  according  to  equation  25: 

a  =  ~  (25) 

i 
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Th  e  association  constant,  K  ,  can  be  derived  from  the  mass- 
balance  equation  utilizing  the  degree  of  dissociation : 

Ka  =  4-=-!  <26> 

a  cy± 

In  equation  26,  y+   is  the  mean  activity  coefficient  on  the 
molar  scale.   For  low  ionic  strength  media  such  as  were 
studied,  the  activity  coefficient  can  be  calculated  with 
the  Debye-Hiickel  law: 

2Az2/T 
f±2  =  10"  1  +  Bavi  (2?) 

where  A  and  B  are  the  usual  Debye-Hiickel  parameters: 

A  =  L8246  x  106 


(28) 
8 


(eT)1-5 

50.29  x  10 

a   _  j 

(eT)2 
This  yields  the  activity  coefficient  on  the  mole  fraction 
scale.   To  obtain  the  coefficient  on  the  molar  scale 
requires  a  conversion  factor: 

y+2  =  ( o e )2  f  2       (29) 

p  +  10   c2M   -.      -  M   ., 

solvent     solute 

where : 

M  =  molecular  weight 

p  =  density  of  solvent 

The  computer  program,  called  "PITTSEQN,"  designed  to 

perform  the  calculations  associated  with  the  Pitts  theory 

is  described  in  more  detail  in  a  following  chapter. 


CHAPTER  III 
EXPERIMENTAL 

Materials 
Solvents 

Sulfolane  was  obtained  from  Eastman,  a  division  of 
the  Eastman  Kodak  Company.   The  material  obtained  was  of 

a  practical  grade   and  was  purified  by  a  double  distillation 

52  54 
at  reduced  pressure.   '     The  first  distillation  was  at 

2  torr  pressure,  stripping  the  boiling  material  with 
nitrogen.   The  second  distillation  was  at  0.5  torr,  from 
solid  potassium  hydroxide.   The  middle  75%  was  collected  in 
each  case,  and  in  the  final  distillation  it  came  off  at 
104°C.   The  product  had  a  specific  conductance  of 
4.09  x  10~8  S  cm-1. 

l-Methyl-2-pyrrolidinone  was  obtained  from  Eastman, 
and  was  reagent  grade.   It  was  purified  by  double  distil- 
lation at  reduced  pressure,  both  distillations  being 
conducted  at  0.5  torr  pressure.   The  second  distillation  was 
from  solid  sodium  hydroxide.   The  middle  75%  was  collected 

in  each  case,  and  the  final  product  had  a  specific 

-7      -1 
conductance  of  3.64  x  10    S  cm 

2-Methoxyethanol  was  obtained  from  Eastman,  and  was 

reagent  grade.   It  was  purified  by  a  single  distillation  at 
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atmospheric  pressure.   The  middle  75%  was  collected. 

When  the  80  weight-%  solvent  was  prepared,  it  was  found 

-7     -1 
to  have  a  specific  conductance  of  1.11  x  10    S  cm 

Conductivity  water  was  prepared  by  the  distillation  of 

water  which  had  previously  been  treated  by  reverse  osmosis 

followed  by  passage  through  a  deionizing  system.   The 

specific  conductance  of  this  water  was  found  to  be  in  the 

-7     -1 
range  on  7-9  x  10    S  cm   ,  quite  reasonable  since  no 

effort  was  made  to  protect  it  from  C0?  in  the  atmosphere. 


Reagents 

Potassium  chloride  (for  determination  of  the  cell 
constants)  was  Mallinckrodt  Analytical  Reagent  Grade, 
purified  by  double  recrystallization  from  conductivity 
water.   Tris  base  ( tris(hydroxymethyl) aminomethane)  was 
obtained  from  Sigma  Chemical  Company  and  was  their  Trizma 
Reagent  Grade.   It  was  purified  by  double  recrystallization 
from  70%  methanol-water .   Sodium  acetate  and  sodium  chloride 
were  Mallinckrodt  Analytical  Reagent  Grade  materials;  they 
were  recrystallized  from  water.   After  recrystallization, 
all  salts  were  dried  -in    vacuo    at  75°C  overnight  and  were 
stored  in  a  desiccator  which  was  charged  with  indicating 
Drierite . 

The  glacial  acetic  acid  used  was  Mallinckrodt  Analytical 
Reagent  Grade  and  was  purified  by  several  fractional 
freezings.   The  final  fraction  was  stored  in  a  desiccator 
charged  with  indicating;  Drierite. 
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Solutions  of  hydrogen  chloride  were  prepared  by  the 
double  distillation,  at  atmospheric  pressure,  of  analytical 
reagent  grade  hydrochloric  acid,  to  obtain  the  constant 
boiling  acid-water  azeotrope.   The  final  concentration  of 
this  azeotrope  was  determined  by  gravimetric  determination 

of  chloride  with  silver  nitrate;  its  composition  was  found 

_3 
to  be  5.27335  x  10    mol  HC1  per  gram  of  solution. 

Hydrogen  chloride  gas  was  prepared  in  two  manners. 
In  the  early  portions  of  the  study,  it  was  prepared  by 
dropping  redistilled  hydrochloric  acid  onto  reagent  grade 
phosphorus  pentoxide.   The  resulting  hydrogen  chloride  gas 
was  forced,  by  a  stream  of  dry  nitrogen,  through  a  tower 
of  indicating  silica  gel,  and  into  a  flask  containing  the 
solvent.   In  the  latter  portions  of  the  study,  the  gas  was 
prepared  by  dropping  reagent  grade  sulfuric  acid  onto 
analytical  reagent  grade  sodium  chloride.   The  resulting 
gas  was  passed  into  the  solvent  via  the  method  described 
above. 

The  concentrations  of  the  acid  solutions  thus  prepared 
were  determined  as  follows.   A  solution  of  analytical 
reagent  grade  sodium  hydroxide  in  water  was  prepared  and  was 
titrated  pH-metrically  against  NBS  primary  standard 
potassium  acid  phthalate.   This  sodium  hydroxide  solution 
was  then  used  in  the  pH  titration  of  the  acid  solution. 

Approximately  0.5  g  of  acid  solution  was  weighed  into 
a  beaker  containing  about  100  ml  of  conductivity  water. 
After  mixing,  this  was  titrated  with  the  sodium  hydroxide 
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solution,  using  a  microburet.   In  the  titrations, 
consistently  sharp  endpoints  were  noted. 

The  concentration  of  the  solution  of  hydrogen  chloride 

-4 
in  sulfolane  thus  determined  was  1.1155  x  10  '  mol  HC1 

g  soln   .   For  the  solution  of  hydrogen  chloride  in  1- 

-4 
methyl-2-pyrrolidinone ,  it  was  found  to  be  5.8359  x  10 

mol  HC1  g  soln-  .   Upon  adding  hydrogen  chloride  gas  to 

sulfolane  a  color  change  was  noticed.   This  intensified 

with  time,  until  a  deep  golden  hue  was  attained.   This 

observation  raised  questions  about  the  purity  of  the 

sulfolane,  and  tests  were  carried  out  to  determine  its 

purity. 

At  first,  infrared  spectroscopy  was  used.   Figure  4 

is  the  resulting  spectrogram.   Superimposed  on  it  is  a 

55 
reference  spectrogram.     Both  the  sulfolane  solvent  and  the 

sulfolane  containing  HC1  gave  identical  spectrograms,  with 

no  evidence  of  a  major  impurity. 

13 
A  second  test  was    C  nuclear  magnetic  resonance 

spectroscopy.   Figure  5  is  the  spectrogram  obtained,  with 

55 
the  reference  spectrogram  superimposed.     Again,  there  was 

no  noticeable  trace  of  other  material. 

Finally,  liquid  chromatography  was  utilized  in  an 

attempt  to  determine  any  trace  impurities.   Figure  6a  shows 

the  chromatogram  of  the  purified  sulfolane;  Figure  6b  shows 

the  chromatogram  of  sulfolane  containing  hydrogen  chloride. 

Neither  shows  any  impurity  in  sufficient  quantity  to  affect 
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the  results  of  this  study.   It  was  concluded  that  despite 
the  formation  of  the  color,  the  sulfolane  utilized  for  this 
study  was  pure. 

Acid  solutions  in  mixed  dipolar  aprotic-water  solvents 
were  prepared  by  one  of  two  methods.   For  those  solvents 
with  a  large  mole  fraction  of  the  nonaqueous  component, 
the  reagent  was  prepared  by  dilution  of  the  stock  solution 
of  acid  in  the  pure  solvent .   Water  was  added  to  bring  the 
solution  to  the  proper  solvent  composition,  and  the 
resulting  solution  was  diluted  to  produce  the  working  stock 
solution. 

For  those  solvents  relatively  rich  in  water,  the 
acid  stock  solution  was  prepared  by  adding  the  nonaqueous 
component  to  redistilled  hydrochloric  acid  to  achieve  the 
proper  solvent  mix.   This  was  then  diluted  to  produce  the 
working  stock  solution. 

All  stock  solutions  for  the  conductance  measurements 
were  prepared  on  a  strictly  weight  basis. 
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Figure  6a 

Chromatogram  of  pure  sulfolane 
Reverse  phase  column 


Figure  6b 

Chromatogram  of  sulf olane-HCl 
Reverse  phase  column 
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Eluant 50/50   CHgCN/II^O 
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Equipment 
The  major  instrument  used  in  this  study  was  a  conduc- 
tance bridge,  in  general  of  a  design  suggested  by  Janz  and 

32 

Mclntyre.     A  schematic  diagram  of  the  bridge  appears  in 

Figure  7.   At  the  heart  of  the  bridge  is  a  General  Radio 
Type  1654-9002  impedance  comparator.   Originally  designed 
to  measure  impedance  at  0.1,  1.0,  10.0  and  100.0  kHz, 
the  instrument  was  modified  by  General  Radio  to  measure  at 
1.0,  2.5,  5.0,  and  10.0  kHz.   The  instrument  utilizes  a 
transformer  bridge  measurement  circuit ,  with  the  trans- 
former bridge  of  toroidal  design.   The  secondary  winding  of 
this  transformer  is  center-tapped,  and  consists  of  two 
identical  wires,  twisted  together  and  wound  on  the  torus. 

The  two  halves  of  this  secondary  act  as  the  ratio  arms  of 

6 
the  bridge   and  are  equal  to  within  one  part  in  10  .   The 

other  two  arms  of  the  bridge  consist  of  the  unknown 

impedance  (the  conductivity  cell)   and  the  standards. 

The  standards  used  in  this  instrument  are  a  General 

Radio  Type  1433-F  decade  resistor  and  a  General  Radio  Type 

1412-BL  decade  capacitor  connected  in  parallel  to  the  decade 

resistor.   The  resistor  has  seven  decades  of  resistance, 

with  steps  of  10  kf2  down  to  steps  of  0.01  ft.   The  capacitor 

has  five  decades  of  capacitance,  ranging  from  steps  of  0.1 

uF  down  to  steps  of  1  pF.   Using  both  of  these  standards,  it 

was  possible  to  measure  both  the  resistance  of  the  test  cell, 

for  conductivity,  and  the  capacitance,  for  dielectric 

constant  measurements.   For  typical  readings,  the  precision 

of  the  standards  is  0.01%. 
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The  null-point  of  the  bridge  was  monitored  with  a 
Keithley  model  191  digital  voltmeter.   This  permitted  far 
more  accurate  zeroing  of  the  meter  (0.001  volt)  than 
visual  observation  of  the  meter  needle. 

To  eliminate  the  problem  of  irreproducible  contact 
resistances  between  the  bridge  circuit  and  the  test  cell, 
and  to  obviate  the  necessity  of  making  corrections  for  lead 
resistances,  the  instrument  was  adapted  to  utilize  a 

four-lead  measurement  system.   This  is  the  same  system 

56a  b 
utilized  with  four-lead  platinum  resistance  thermometers 

and  requires  two  measurements  with  alternate  sets  of  leads. 

The  mean  of  the  two  measurements  effectively  averages  out 

resistances  inherent  to  the  leads  and  the  connections. 

To  further  cut  down  on  problems  associated  with 
connection  of  the  cell  to  the  bridge,  printed  circuit  boards 
and  edge  connectors  were  used.   Standard  equipment  with  15 
contacts  was  used,  allowing  three  contacts  per  lead,  with 
one  contact  between  live  connections  as  a  spacer.   The 
multiplicity  of  contacts  per  lead  insured  good  and 
reproducible  contact. 

The  conductance  cells  used  in  this  study  were  of  the 
type  designed  by  Daggett,  Bair,  and  Kraus,    and  were 
obtained  from  Beckman  Instruments.   The  constants  determined 
for  the  cells  were  1.05393  and  0.120695  era"  .   The  cell 

appears  in  Figure  8.   The  leads  are  held  well  separated  to 

24  ^7 
help  reduce  the  Parker  effect.   '     This  effect  is  caused 

by  the  electrical  leads  passing  through  the  electrolyte 


Figure  8 
Conductance  cell 
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solution,  and  gives  rise  to  a  capacitative  effect  causing 
a  frequency  dependence  in  the  measured  resistance. 

Another  source  of  error  in  conductance  measurements 

related  to  the  measuring  cell  is  the  Soret  effect.   This 

33 
has  been  discussed  in  depth  by  Stokes.     Briefly,  this 

effect  is  due  to  a  slow  thermal  diffusion  of  the  electrolyte 
within  the  electrode  compartment.   The  design  of  the  flask- 
type  dilution  cell  used  in  this  study,  coupled  with  constant 
stirring,  minimizes  this  effect. 

The  electrodes  of  the  cell  were  shiny  platinum,  with 
a  circular  cross  section.   They  were  spot  welded  to 
platinum  wire,  which  in  turn  was  welded  to  tungsten  wire  to 
achieve  a  seal  with  the  glass.   Heavy  gauge  copper  wire  was 
welded  to  the  tungsten  to  form  the  cell  leads,  and  the 
copper  wire  was  doubled  at  the  top  of  the  cell  to  produce 
the  four-lead  measurement  system. 

With  two  cells  in  the  temperature  bath  simultaneously, 
held  in  close  proximity  to  one  another,  any  shaking  could 

result  in  mutual  damage.   For  this  reason,  a  special  filler 

57 
cap  was  designed    to  facilitate  additions  to  the  cell 

without  shaking.   The  filler  cap  and  its  companion  weight 

buret  appear  in  Figure  9.   By  adjusting  the  position  of  the 

3-way  stopcock,  additions  could  be  made,  and  then  solution 

forced  up  the  fill-tube  into  the  upper  chamber.   This 

permitted  the  washing  down  of  any  splashes  made  during  the 

addition  and  facilitated  solution  mixing.   Both  the  weight 

buret  and  the  filler  cap  were  manufactured  by  the 


Figure  9 
Filler  cap  and  weight  buret 
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departmental  glassblower  and  utilize  standard  ground-glass 
joint  connections. 

The  constant-temperature  bath  utilized  in  this  study 
was  constructed  by  the  departmental  machine  shop  out  of 
stainless  steel,  with  all  seams  welded.   It  was  insulated 
with  styrofoam  and  was  supported  by  an  aluminum  box  frame 
mounted  on  a  welded  angle-iron  cart.   The  bath  was  filled 

with  BP  food-grade  white  oil,  which  further  helped  to 

95 
minimize  the  Parker  effect  mentioned  previously."    The  oil 

was  circulated  by  two  American  Instrument  Company  1.3  A 

circulating  pumps,  to  assure  temperature  uniformity 

throughout  the  bath. 

Temperature  control  was  achieved  with  a  Yellow 
Springs  Instrument  Company  model  72  proportional  temperature 
controller.   The  usual  three-quarter  turn  potentiometer 
used  for  regulation  of  f ractions-of-a-degree  was  replaced 
with  a  ten-turn  potentiometer  to  provide  far  more  precise 
control  of  the  temperature  setting.   The  controller  was 
used  with  a  YSI  thermistor  probe  and  was  used  to  drive  two 
heaters.   One  was  a  standard  500-watt  Vycor  sheathed 
immersion  heater.   The  other  was  constructed  of  nichrome 
heating  wire  wound  around  a  plexiglas  frame;  it  was  placed 
to  take  full  advantage  of  the  flow  from  the  two  circulating 
pumps. 

Immersed  in  the  bath  were  a  number  of  coils  of  copper 
tubing.   Through  these  coils  was  circulated  water  from  an 
external  bath.   This  bath  was  cooled  by  a  Blue-M 
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Constant-Flow  model  PCC-24A-2  cooling  unit ,  which  has 
thermostatic  control.   This  external  water  bath  was  kept 
approximately  3°C  cooler  than  the  temperature  of  the  main 
bath;  as  a  result,  the  heaters  in  the  main  bath  were  never 
completely  shut  off.   This  allowed  constant  regulation  of 
the  temperature  of  the  bath  rather  than  the  cycling  usually 
associated  with  proportional  temperature  controllers.   The 
water  from  the  external  bath  was  circulated  by  an  American 
Instrument  Company  1.3  A  circulating  pump,  and  the 
regulation  of  temperature  in  the  main  bath  by  the  apparatus 
described  above  was  ±0.002°C. 

Temperatures  were  monitored  with  a  Hewlett-Packard 
model  2801A  digital  quartz  thermometer  and  quartz  probe. 
This  digital  thermometer  was  calibrated  across  the  entire 
temperature  range  of  interest,  and  the  appropriate  correction 
factors  were  determined.   It  was  calibrated  against  a  Leeds 
and  Northrup  platinum  resistance  thermometer,  model  #8160. 
The  resistance  of  the  resistance  thermometer  was  measured 
with  a  Leeds  and  Northrup  type  G-2  Muller  bridge,  model  8069, 
which  allowed  measurement  of  thermometer  resistance  to 
0.0001  ft.   A  response  printout  from  Leeds  and  Northrup 
permitted  the  determination  of  temperature  measured  directly 
from  the  resistance  reading. 

To  further  cut  down  on  heat  loss  from  the  bath,  and 
to  prevent  debris  from  entering  the  oil,  a  plexiglas  lid 
for  the  bath  was  constructed  by  the  departmental  machine 
shop.   This  lid  had  a  hinged  door  over  the  cell  area  to 
allow  access  to  the  cells  therein. 
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The  cells  were  stirred  constantly  by  a  Poly  Science 
Corporation  model  RZR-10  stirrer  motor,  which  permitted 
variable  stirrer  speed.   The  motor  was  connected  via  a 
driveshaft  and  two  Sears  Craftsman  right-angle  drives  to 
a  Pic  chain-and-sprocket  drive.   This  drive  turned  two 
large  teflon  stirbars  positioned  under  the  two  cells.   The 
temperature  bath  and  its  associated  components  are  diagrammed 
in  Figure  10. 

Dielectric  constants  were  determined  using  a  Balsbaugh 
Laboratories  model  100-T3  dielectric  constant  cell.   This 
cell  had  a  three-terminal  arrangement  and  was  constructed 
of  cylindrical  nickel-steel  electrodes.   A  special  wiring 
harness  was  constructed  to  allow  connection  to  the  four-lead 
measuring  system  of  the  conductance  bridge.   The  cell  was 
thermostated  in  the  constant  temperature  bath,  and  after 
thermal  equilibrium  was  reached  no  appreciable  drift  in 
capacitance  was  noted. 

Viscosities  were  determined  with  a  Canon-Fenske  type 
viscometer,  calibrated  with  conductivity  water  at  the 
various  temperatures  studied. 

Densities  were  determined  with  three  pycnometers  of 
the  Gay-Lussac  type,  obtained  from  the  Fisher  Scientific 
Company.   The  pycnometers  were  calibrated  with  conductivity 
water  at  the  temperatures  studied,  and  the  three  readings 
were  averaged  together  to  obtain  the  mean  density. 
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CHAPTER  IV 
DATA  ACQUISITION  AND  TREATMENT 

Experimental  Technique 
A  standard  procedure  for  conductance  measurements  was 
adopted  and  adhered  to  for  this  study.   An  amount  of 
solvent  (usually  about  350  g)  was  weighed  into  the 
conductance  cell,  which  was  then  placed  into  the  constant 
temperature  bath.   After  thermal  equilibrium  was  reached 
(in  about  an  hour),  the  resistance  of  the  solvent  was 

measured.   Since  solvent  resistances  generally  range  from 

5      6 
10   to  10  Q,    and  since  the  decade  resistance  standard  had 

a  maximum  resistance  of  99,999.99  ft,  it  was  found  necessary 

to  employ  a  resistive  shunt.   This  shunt,  in  parallel  to 

the  conductance  cell,  was  a  precision  resistor  whose 

resistance  was  evaluated  on  the  bridge  as  10,006.8  ft. 

Response  was  found  to  be  flat  across  the  frequency  range 

of  the  measurements.   This  allowed  measurement  of  the  very 

high  resistances  according  to  the  equation: 

R  R 

Ra  =  R~ r-H-  (30) 

s     m 

where  R   is  the  solution  resistance,  R   is  the  measured 
a  '   m 

resistance,  and  R   is  the  shunt  resistance. 

s 


-46- 


-47- 

Once  the  solvent  resistance  was  determined,  the  cell, 
with  contents  continuously  stirred,  was  allowed  to  remain 
undisturbed  and  the  resistance  was  monitored  as  a  function 
of  time.   If  the  solvent  were  to  leach  any  impurity  out  of 
the  glass  of  the  cell,  the  resistance  would  decrease  with 
time.   This  was  not  noted  for  any  of  the  solvent  systems  in 
this  study. 

A  concentrated  solution  of  the  electrolyte  under  study 
was  prepared  in  the  solvent  of  interest.   Weighed  amounts 
of  this  stock  solution  would  be  added  incrementally  to  the 
cell  via  the  filler  cap.   After  each  addition,  the  fill 
tube  would  be  rinsed  with  the  solution  in  the  cell  by 
forcing  the  solution  up  into  the  fill  tube  with  the  rubber 
bulb.   After  equilibrium  was  reached,  as  indicated  by 
stability  of  the  measured  resistance  with  time,  the  solution 
resistance  was  recorded  and  another  addition  of  stock 
solution  was  made.   Usually,  12  such  solutions  were 
measured,  and  after  the  last  measurement  the  density  of  the 
final  solution  was  measured. 

Resistances  were  measured  with  a  test  voltage  of  0.3  V 
at  three  frequencies:   2.5,  5.0,  and  10.0  kHz.   These 
resistances  were  then  extrapolated  vs.  a  function  of 
frequency  to  eliminate  the  frequency  dependence,  as  described 
in  a  following  section. 

Solution  equilibrium,  for  the  2-methoxyethanol/water 
solvent  system,  was  found  to  be  attained  in  about  15 
minutes.   For  the  sulf olane/water  system,  equilibrium  took 
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up  to  one  hour,  probably  due  to  the  very  slow  solvation  of 
the  cation  by  the  solvent.   For  the  l-methyl-2-pyrrolidinone/ 
water  system,  equilibrium  was  achieved  in  about  15  minutes. 
Quickest  of  all  was  the  system  of  Tris-HCl  in  water,  with 
equilibrium  in  less  than  5  minutes. 

Making  the  conductance  measurements  in  the  manner 
described  above  gave  rise  to  a  set  of  concentration  (in  moles 
of  solute  per  gram  of  solution)  vs.  resistance  points. 
Treatment  of  these  data  is  described  in  a  following  section. 

Dielectric  constant  measurements  were  greatly 
facilitated  by  the  use  of  a  three-terminal  dielectric 
constant  cell.   The  bridge  was  balanced  (shunt  in  place) 
with  the  cell  disconnected.   With  the  cell  connected,  the 
capacitance  of  the  empty  cell  (i.e.,  with  air  as  the 
dielectric  medium)  was  measured.   The  cell  holder  was  then 
filled  with  the  solvent  of  interest  and  placed  into  the 
constant  temperature  bath.   The  cell  was  inserted,  and  the 
solution  was  allowed  to  come  to  thermal  equilibrium.   The 

capacitance  was  then  measured,  and  the  dielectric  constant 

58 
was  very  simply  calculated  with  the  equation: 

C 
e  =  ^  (31) 

a 

where  C   is  the  solvent  capacitance  and  C   is  the 
s  ^  a 

capacitance  of  the  empty  cell. 

Cell  constants  were  determined  utilizing  the  method 

59 
described  by  Lind,  Zwolenik,  and  Fuoss.     Utilizing  data 

available  for  the  conductance  of  aqueous  potassium  chloride, 
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they  derived  an  equation  to  give  the  equivalent  conduc- 
tance, A,  for  a  given  concentration  of  KC1 : 

A  =  149.93  -  84.65/c  +  58.74c  log  c  +  198.4c        (32) 
where  c  is  the  molar  concentration  of  KC1 .   They  found 

this  equation  to  give  excellent  results  when  compared  to  the 

26 
Jones  and  Bradshaw    conductance  standard,  probably  the 

most  widely  accepted  standard  for  conductance  measurements. 

The  method  was  utilized  by  making  a  conductance  run 
using  water  as  the  solvent  and  adding  increments  of  a 
concentrated  KC1  stock  solution.   Once  the  concentrations 
were  converted  to  molarity  and  the  resistances  extrapolated 
to  remove  frequency  dependence,  the  cell  constants  were 
calculated. 

As  was  shown  in  a  previous  chapter,  the  resistance  of 
a  solution  is  directly  relatable  to  the  specific 
conductance : 

k  =  ~  (33) 

s 

The  solution  resistance  is  the  measured  resistance,  after 

extrapolation,  corrected  for  the 'resistance  of  the  solvents: 

J_  =  J_ L_  (34) 

R       R       R  ^     ' 

rt      rt      rt    n 

s    m    solv 
The  specific  conductance  may  also  be  related  to  the 
equivalent  conductance: 

K  =  lOOO  (35) 
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Combining  equations  33  and  35  we  obtain  an  expression 
relating  cell  constant,  solution  resistance  and  equivalent 
conductance : 

A   -   10009  flR^ 

A  -  -ir5-  (36) 

s 

where  c  is  the  concentration  in  mol  L 

By  plotting  A,  calculated  with  equation  32,  vs.  the 
quantity  (R  c)   ,  a  line  is  obtained  whose  slope  is  10000. 

In  this  study,  two  separate  calibration  runs  were 
made,  with  a  total  of  10  measurements  for  each  of  the  two 
conductance  cells.  The  cell  constants  were  determined  by 
linear  regression  analysis  to  be  1.05393  cm-  and  0.120695 
cm  .  In  both  -cases ,  the  correlation  coefficient  for  the 
regression  was  greater  than  0.99999,  indicating  extremely 
good  reproducibility  between  runs. 

Recently,  an  alternate  equation  for  the  determination 

of  the  equivalent  conductance  of  potassium  chloride  solutions 

60 
has  been  proposed  by  Justice.     His  equation  has  the  form: 

A  =  149.87  -  94.87/c  +  58.63c  log  c  +  224.1c  -  254.5c1,0  (37) 
It  can  be  seen  that  this  equation  has  a  similar  form  to 

the  equation  of  Fuoss  et  al .  but  is  extended  to  include 

1  5 
the  c  '   term.   By  doing  so,  the  applicability  of  the 

equation  is  extended  up  to  approximately  0.04  M  KC1,  as 

opposed  to  0.012  M  KC1  for  the  Fuoss  equation. 

44 
A  more  recent  paper  by  Fuoss  and  Hsia    contains  a 

modification  of  the  original  equation  which  incorporates 

the  c     term: 
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A  =  149.936  -  94.88/c  +  25.48c  In  c  +  221.0c  -  299c1-5   (38) 
As  can  be  seen,  the  constants  differ  slightly  from  the 
equation  proposed  by  Justice.   It  is  stated  that  the  new 
Fuoss  equation  will  provide  agreement  with  the  Jones  and 
Bradshaw  0.1  demal  KC1  standard,  thus  extending  the  concen- 
tration range  over  which  the  cell  can  be  calibrated. 

Since  this  study,  however,  was  done  at  concentrations 
on  the  order  of  0.01  M  maximum,  the  original  Fuoss  equation 
was  used  for  the  determination  of  the  cell  constants.   At 
such  low  concentrations,  the  difference  between  the  A 
values  generated  by  that  equation  and  the  values  generated 
by  the  two  later  equations  is  less  than  one-half  of  one 
percent . 
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Preliminary  Data  Treatment 

As  has  been  mentioned  previously,  concentrations  were 
determined  by  weight,  and  solution  additions  were  by 
weight.   This  resulted  in  experimental  concentrations  in 
units  of  mol  g  soln   .   Since  equivalent  conductance 
is  in  terms  of  molarity,  it  was  necessary  to  convert  to 
molar  concentrations.   This  was  done  through  use  of  the 
density  of  the  solution. 

Both  the  solvent  density  and  the  density  of  the  final 
most  concentrated  conductance  solution  were  determined.   The 
densities  of  the  more  dilute  solutions  were  assumed  to  lie 
on  a  straight  line  connecting  the  two  experimental  points. 
The  concentration  in  mol  g  soln    and  the  molar  concentration 
may  be  related  by  the  equation: 

c(M)  =  1000m(pslm  +  Psoln)  (39) 

where  m  is  the  concentration  in  mol  g  soln   ,  p  ,  is  the 

slope  of  the  line  connecting  the  two  experimental  densities, 

and  p   ,   is  the  solvent  density.   The  quantity  p  -,  m  +  P  -, 
solv  J  i  J    'si      solv 

is  in  fact  the  density  of  the  solution  whose  concentration 
is  m. 

Resistances  were  measured  at  2.5,  5.0,  and  10.0  kHz 
and  were  then  extrapolated  to  infinite  frequency  by  use  of 

the  method  endorsed  by  Hoover  for  shiny  platinum  electrodes 

3^ 
such  as  were  used  in  this  study.  *■   The  extrapolation  is 

based  on  an  equation  of  the  form: 

R2 

R  ,    =   R   .      .       + ft  (40) 

measured  ohmic         .2  v       ' 

1    +    auj 
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where  : 


F  -  1 

2    _   2 
L0o   -  Foj„ 


and : 


2      2  D 
*    "l   "  C°3  -  2  "  Rl 

P   =  .... . ~» 

2      2  R0  —  R-, 

03-.    -  OJp     3      1 

with  co  =  2ttv.   At  infinite  frequency,  the  term  on  the  right- 
hand  side  of  equation  40  drops  out,  leaving  the  true  ohmic 

resistance.   The  measured  resistance  is  plotted  vs. 

2  -1 
( 1  +  ato  )   ;  the  intercept  is  the  ohmic  resistance. 

Usually,  the  measured  resistance  was  found  to  decrease 

with  increasing  frequency,  which  is  considered  normal. 

However,  when  dealing  with  very  high  resistances,  such  as 

those  associated  with  pure  solvents,  a  reversal  of  this 

trend  was  noted;  the  resistance  increased  with  frequency. 

61 
This  phenomenon  has  been  described  by  Mysels  et  al .   .  and 

is  attributed  to  a  leakage  to  ground.   Their  recommendation 

is  that  the  resistance  be  extrapolated  to  zero  frequency  vs. 

2 

co  .   However,  their  recommendations  are  for  frequencies 

of  measurement  less  than  0.5  kHz,  far  lower  than  those 
employed  in  this  study. 

Where  the  resistances  were  found  to  be  sufficiently 
high  to  cause  anomalous  frequency  dependence  (pure  solvents, 
HC1  in  sulfolane,  acetic  acid  in  2-methoxyethanol/water) 
the  resistance  values  at  10  kHz  were  used  for  further 
calculation  of  conductance  parameters. 
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The  preliminary  data  treatment  to  this  point  was 
performed  with  a  single  computer  program,  called  "FREQEXTR," 
written  in  the  Level  II  Basic  language  to  utilize  a  Radio 
Shack  TRS-80  microcomputer,  operating  under  the  NewDos+ 
disk  operating  system,  with  32  KBytes  of  random  access 
memory  (RAM).   A  copy  of  the  program  appears  in  the  Appendix. 
After  inputting  the  various  parameters  (density  slope, 
solvent  density,  cell  constant)  the  resistances  and  the 
corresponding  frequencies  were  input,  along  with  the 
experimental  concentration.   The  computer  supplied: 

1)  The  resistance  at  infinite  frequency; 

2)  The  concentration  in  mol  L   ;  and 

3)  The  equivalent  conductance,  calculated  via  equation  36. 
Solvent  corrections  were  always  applied. 

Treatment  of  the  experimental  data  with  this  program 
yielded  a  set  of  molar  concentration  vs.  equivalent 
conductance  data  points  which  were  then  processed  as  des- 
cribed below. 


-55- 

Data  Processing 
The  conductance  vs.  concentration  data  were  processed 
with  a  computer  program  designed  to  utilize  the  full  Pitts 
equation  described  previously.   The  program  was  initially 
written  in  the  Fortran  IV  language  to  utilize  the 

facilities  of  the  Northeast  Regional  Data  Center  on  the 

62 

campus  of  the  University  of  Florida.  "   With  the  acquisition 

of  the  microcomputer  system  described  previously,  the 
program  was  translated  into  the  Basic  language;  at  the  same 
time,  the  program  was  streamlined  and  rendered  more 
efficient.   The  Fortran  IV  subroutine  "EXPI"  for  the 
evaluation  of  exponential  integrals  was  translated  in  full. 
All  calculations  were  carried  out  in  double  precision. 

The  program  is  designed  to  utilize  an  iterative 

63  64 
procedure  devised  by  Ives,   '    using  the  following 

expression  for  the  dissociation  constant,  K, : 

2   2 

ci  y  + 

K,  =  — —  (41) 

d     c  v  ' 

u 

where  the  subscripts  i  and  u  refer  to  ionic  and  nondis- 

sociated  quantities,  respectively.   Since  the  degree  of 

dissociation,  a,  may  be  expressed  as: 

a=A  =  fA  (42) 

l 

both  c.  and  c   mav  be  expressed  in  terms  of  total  concen- 
1       u  v 

tration  of  electrolyte  and  the  degree  of  dissociation: 


(43) 
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c .  =  ac 

1 

c   =  (1  -  a)c 

u 

Combining  the  preceding  four  equations,  one  arrives  at 
an  expression  for  K,  in  terms  of  A: 

A2cn2 
Kd  "   A. (A."-   A)  <44> 

11  ' 

In  general,  A.,  the  ionic  equivalent  molar  conductance, 
may  be  expressed  as  a  function  of  A  ,  the  limiting 
equivalent  conductance: 

A.  =  A°  +  A  (45) 

l         r 

where  A   is  an  expression  for  altering  A   to  obtain  A. ; 

A  will  vary  with  the  theory  being  applied. 

K  ,  may  now  be  expressed  as : 

.2    2 

A  cy  + 

K   = (46) 

a    A.(A°  +  A  )  -  A 
iv       rJ 

or: 

2c   2 

A  -  Ar  =  A°  -  A_J^L_  (47) 

i  d 

If  we  now  set  up  two  variables,  X  and  Y,  such  that: 

Y  =  A  -  A 
r 

A2cy,2 
X  =  -A7=- 

1 

we  may  reduce  the  equation  to  the  form: 


(48) 


Y  =  A°  -  £-  (49) 

Kd 
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With  both  X  and  Y  known,  this  reduces  to  a  simple  linear 
least-squares  procedure  for  obtaining  A   (the  intercept), 
and  K,  (-slope)   .   Since  X  and  Y  depend  on  A  ,  Ives 
proposed  an  iterative  procedure  whereby  an  estimate  for 
A   is  provided.   A  set  of  X  and  Y  values  is  compiled  using 
this  estimate,  a  least-squares  calculation  is  done,  and  a 
new  value  for  A   is  obtained.   This  is  then  used  to  generate 
new  values  for  X  and  Y,  and  so  forth,  until  no  change  in  A 
is  noted.   This  usually  takes  less  than  five  iterative 
steps . 

In  his  original  treatment ,  Ives  used  the  Onsager 
limiting  law  to  obtain  A. ;  an  exact  solution  for  this 
quantity  can  thus  be  obtained.   When  using  the  Pitts  theory 
of  conductance,  however,  this  is  not  the  case.   From 
Equation  45: 

A.  =  A°  +  A  (45) 

l         r  v  ' 

In  the  Pitts  equation,  A   is  a  transcendental  function  of 
m       .   r 

A.  and  a  direct  solution  is  not  possible.   Hence,  the 
program  was  designed  to  iterate  to  obtain  the  A   value. 
The  program  is  designed  to  take  the  estimate  of  A 
as  the  first  approximation  of  A . .   It  then  calculates  the 
various  parameters  needed  (a,  c . ,  Pitts  terms,  etc.)  and 
using  these  calculates  the  Pitts  A   term.   Using  this,  it 
calculates  A.  and  compares  it  to  the  original  estimate.   If 
the  difference  is  greater  than  0.001,  it  takes  the  new  A. 
value  and  goes  through  the  iteration  process  again.   This 
is  continued  until  convergence  occurs. 
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The  iterative  process  described  above  takes  place  for 

each  data  point;  once  it  is  concluded,  the  X  and  Y  values 

are  generated,  and  the  program  proceeds  to  the  next  point. 

When  all  points  have  been  processed,  the  program  performs 

the  least-squares  calculation  to  obtain  the  values  for 

A   and  K , .   In  this  study,  a  slight  alteration  was  made  to 

have  the  program  calculate  the  association  constant 

K   (=  77—),    rather  than  the  dissociation  constant, 
a    nd 

This  procedure  was  carried  out  for  five  iterations; 

values  of  A   and  K   have  usually  converged  by  the  fourth. 

a 

The  complex  point-by-point  iterations  repeated  five  times 
explain  why  the  microcomputer  execution  time  for  this 
program  runs  in  excess  of  3.5  hours  for  a  typical  set  of 
conductance  data!   The  performance  of  the  program  has  been 

tested  and  found  to  yield  results  in  very  good  agreement 

57 
with  those  considered  standard. 

As  mentioned  in  a  previous  chapter,  the  full  Pitts 

equation  has  been  expanded  by  Fernandez-Prini  and  Prue  to 

give  an  equation  of  the  familiar  form: 


A 


1.5 


=  A°  -  S/c  +  Ec  log  c  +  J  c  -  J2  (23) 


During  his  visit  to  the  laboratory,  Dr.  Fernandez-Prini 
graciously  helped  to  prepare  a  computer  program  designed  to 
utilize  this  form  of  the  conductance  equation,  for  both  the 
Pitts  and  Fuoss-Hsia  theories.   This  program  was  also 
adapted  for  use  with  associated  electrolytes  by  replacing  c 
with  ac.   Since  the  conductance  due  to  free  ions  (the  ionic 
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conductance)  at  concentration  ac  is  given  by  the  following 
relation : 

A  =  Ai(ac)  -  Ay±2Kaac  (50) 

derived  from  equation  26,  the  final  form  of  the  conductance 
equation  for  associated  electrolytes,  combining  equations 
23  and  50  becomes : 

A  =  A°  -  S/oc  +  Eac  log  ac  +  J^ac  -  J2(ac)1,0  -  K  Ay+2ac  (51) 

As  mentioned  earlier,  the  S  and  E  terms  are  identical  for 
both  theories;  the  J-.  and  J?  terms  differ.   The  program 
using  this  form  of  equation  has  several  advantages  over  the 
program  utilizing  the  full  Pitts  equation.   It  can  give  us 

the  results  using  both  theories;  using  the  least-squares 

65 
iterative  method  of  Kay  '  it  operates  in  single  precision 

and  takes  only  a  few  minutes  to  complete  the  calculations. 

However,  it  has  one  major  disadvantage. 

In  the  course  of  the  study,  it  was  found  that  the 

program  utilizing  the  expanded  equations  broke  down  for  the 

2-methoxyethanol/water  and  l-methyl-2-pyrrolidinone  systems. 

The  similarity  here  is  that  both  systems  have  dielectric 

constants  of  about  32,  and  consequently  Bjerrum's  distance 

o 

is  about  8.7  A.   This  means  that  k<x  -  0.45  at  a  concentration 

40  60  66 
of  0.01  M,  and,  as  has  been  pointed  out,  in  this 

region  the  expansions  tend  to  break  down.   This  was 

manifested  by  a  divergence  in  A   values  rather  than  the 

expected  convergence. 
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For  this  reason,  among  others,  the  full  Pitts  equation 
program,  "PITTSEQN"  as  described,  was  used  to  analyze  the 
data  in  this  study.   It  might  be  pointed  out  that,  where  the 
expanded  equation  program  worked,  it  was  quite  useful  for 
obtaining  a  close  approximation  of  A   to  use  as  the  estimate 
for  A°  in  "PITTSEQN";  this  cut  down  initial  iteration  time. 
Copies  of  both  computer  programs,  "PITTSEQN"  and  "EXPEQNS," 
in  the  TRS-80  Level  II  Basic  language,  appear  in  the  Appendix, 

In  both  programs  it  will  be  noted  that  when  the 
equivalent  conductance  is  input  it  is  multiplied  by  the 
factor  0.999505.   This  is  to  convert  the  conductance  in 
international  ohms,  as  measured,  into  the  absolute  ohms 
which  are  standard. 

For  those  systems  studied,  where  the  solvent  conduc- 
tance was  an  appreciable  fraction  of  the  total  conductance 
(i.e.,  where  the  association  constant  was  exceedingly  high), 
neither  program  was  able  to  perform  the  extrapolation  to 

infinite  dilution.   Rather,  another  program  was  written  to 

42 
utilize  the  Shedlovsky-Kay  extrapolation  method.     This 

was  devised  for  use  with  data  for  very  weak  acids  in  mixed 
solvent  systems  (e.g.,  acetic  acid  in  methanol-water  mix- 
tures), and  depends  on  obtaining  an  estimate  for  A  ,  which 
cannot  be  directly  obtained  from  the  data.   This  is  done  by 
using  the  Kohlrausch  law  of  independent  migration,  and  will 
be  discussed  in  depth  in  the  following  section.   The  program 
for  this  extrapolation,  "SHEDKAY,"  may  be  found  in  the 
Appendix. 


CHAPTER  V 
RESULTS 

This  study  may  be  divided  into  three  main  sections: 
first,  an  examination  of  the  conductimetric  behavior  of 
HC1  in  sulfolane  and  sulf olane-water  solvents  of  varying 
compositions  at  various  temperatures;  second,  a  study  of 
the  conductimetric  behavior  of  HC1  in  l-methyl-2- 
pyrrolidinone  (NM^Py)  and  NMpPy-water  solvents  of  varying 
compositions;  and  third,  and  examination  of  the  behavior  of 
various  electrolytes  in  the  solvent  consisting  of  80  weight-% 
2-methoxyethanol  and  20  weight-%  water  in  order  to  determine 
the  amount  of  ion  pairing  present  for  the  electrolytes. 

In  order  to  avoid  the  inclusion  of  large  amounts  of 
data  in  this  chapter,  the  actual  conductance  results  have 
been  placed  in  Appendix  A,  which  may  be  found  following 
Chapter  VII  of  this  dissertation. 

Table  I  is  a  summary  of  the  calculated  conductance 

parameters  for  HC1  in  sulfolane  and  sulf olane-water  solvents 

at  25,  30,  and  40°C.   The  solvent  composition  is  expressed 

as  mole  fraction  of  sulfolane  (X  ) ;  A   is  the  limiting 

equivalent  conductance,  K   is  the  experimental  ion  pairing 

a 

constant,  a    is  the  value  used  for  the  distance  of  closest 
approach  parameter  in  the  Pitts  conductivity  theory,  and 
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aA_A   ^s  ^he   standard  deviation  between  the  experimental 
c 

and  calculated  equivalent  conductances. 

The  distance  of  closest  approach  parameter,  a,    was, 
in  this  study,  set  equal  to  the  Bjerrum  critical  distance; 
this  is  discussed  in  the  following  chapter. 

Table  II  presents  the  calculated  conductance  parameters 
for  HC1  in  NM„Py  and  NM2Py-water  solvents.   The  symbols 
have  the  same  meanings  as  in  Table  I. 

Table  III  presents  the  calculated  conductance  results 
for  several  electrolytes  in  the  80  W%  2-methoxyethanol-water 
(80  W%  MetOH-H20)  solvent.   These  electrolytes  were  chosen 
in  view  of  their  general  importance  in  acid-base  chemistry, 

and  to  tie  in  with  other  research  under  way  at  this 

71 
laboratory.     Furthermore,  for  a  comparative  evaluation  of 

the  results  obtained  in  80  W%  MetOH-H20  with  those  in  pure 

water,  the  behavior  of  aqueous  solutions  of  Tris-HCl 

( tris(hydroxymethyl )  aminomethane  hydrochloride)  was 

72 
studied.     The  results  of  this  study  are  summarized  in 

Table  IV. 

In  this  case,  various  trial  values  of  a   were  used  to 

obtain  a  minimum  in  the  standard  deviation  between  the 

experimental  and  calculated  conductances.   The  value  of 

a   which  gave  the  minimum  in  o._.      was  that  selected  for 

c 
use  in  generating  the  conductance  parameters  listed  in 

Table  IV. 
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TABLE  IV 
Calculated  conductance  parameters  for  Tris-HCl  in  water. 

T  A°  K  a  a.  . 

a  A-Ac 

(°C)     (S  cm2  mol-1)     (mol-1  L)     (cm  x  10~8) 
25         106.07  1.01  2.60        0.225 

37         133.65  0.99  3.75        0.098 


CHAPTER  VI 
DISCUSSION 

Solvent  Properties 
General 

In  general,  electrolytic  conductance  will  depend,  to 
a  large  degree,  on  the  properties  of  the  solvent.   If 
the  solvent  is  able  to  solvate  strongly  the  ions  present , 
this  will  aid  the  dissociation  process  and  increase  the 

conductance  of  the  solution. 

39 

There  are  two  basic  types  of  electrolytes.     In  the 

first,  the  ions  are  already  present,  as  in  inorganic  salts, 
and  dissociate  when  dissolved.   In  the  second  type,  ions 
are  formed  by  reaction  between  the  solvent  and  the  solute; 
for  example,  the  reaction  which  takes  place  when  an  organic 
acid  is  dissolved  in  water: 

RCOOH  +  H20  t    RCOO"  +  H30+  (52) 

The  most  obvious  force  between  a  neutral  solvent 
molecule  and  an  ion  is  the  electrostatic  force  arising  from 
the  interaction  of  the  charge  on  the  ion  with  the  permanent 
dipole  moment  of  the  molecule.   If  the  solvent  molecule  has 
no  permanent  dipole  moment ,  one  may  be  induced  by  the 
presence  nearby  of  the  charged  ion.   Of  a  similar  nature  are 
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the  dispersion,  or  London  forces.   These  are  the  forces 
responsible  for  the  cohesion  of  molecules  lacking  charge 
or  dipole  moments  and  are  due  to  random  formation  of  dipoles 
due  to  fluctuations  in  the  electron  clouds  of  the  molecules. 

Aside  from  these  non-specific  forces,  there  are  several 
types  of  specific  forces  aiding  the  solvation  process.   A 
good  example  of  such  is  hydrogen  bonding.   An  example  is 
the  dissolution  of  an  acid  in  water,  as  depicted  in  equation 
52.   It  is  this  hydrogen  bonding,  and  the  solvent  structure 
that  it  gives  rise  to,  that  leads  to  the  proton  jump 
mechanism  which  causes  abnormally  high  conductances  for 
strong  acids  in  water.   This  will  be  discussed  in  a  later 
section . 

Another  of  the  factors  whereby  the  solvent  influences 
the  conductance  of  an  electrolyte  therein  is  the  viscosity. 
Quite  obviously,  as  the  viscosity  of  the  solvent  increases, 
the  drag  on  moving  ions  increases,  and  the  conductance  will 
thus  decrease.   Yet  another  factor  is  the  solvent  dielectric 
constant.   From  Coulomb's  Law,  it  is  seen  that,  as  the 
dielectric  constant  decreases,  the  electrostatic  forces  of 
attraction  increase;  thus,  the  mutual  attraction  between 
oppositely  charged  ions  increases,  increasing  ion  association 
and  decreasing  conductance. 

In  view  of  the  above,  it  is  interesting  to  examine  the 
properties  of  the  solvent  systems  utilized  in  this  study 
and  to  see  if  any  conclusions  can  be  reached  concerning 
solvent  properties  and  conductance. 
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Sulfolane  System 

The  physical  properties  of  sulfolane  and  sulfolane- 
water  solvents  are  summarized  in  Table  V.   They  are 
primarily  experimental  values,  supplemented  where 
appropriate  by  values  from  the  literature. 

TABLE  V 
Physical  properties  of  sulfolane-water  solvents. 


Temp. 
(°C) 

X   ,  . 
sulf 

K 

(S  cm-1) 
(  x  10"7) 

Density 
(g  mL"1) 

Viscosity 
(cP) 

Dielectric 
Constant 

25 

0.85 

0.735 

1.2564 

7.490 

45.23 

0.75 

2.658 

1.2492 

6.772 

46.77 

0.50 

4.5 

1.2290 

4.583 

52.28  a 

0.25 

6.3 

1.1726 

2.738 

59.70  a 

0.00 

7.0 

0.99707 

0.890 

78.36  b 

30 

1.00 

1.260 

1.2623 

10.30 

43.30 

0.85 

0.298 

1.2535 

7.236 

44.69 

0.75 

20.88 

1.2464 

5.994 

46.26 

0.50 

2.740 

1.2214 

3.797 

51.25 

0.25 

3.610 

1.1698 

2.287 

57.57 

0.00 

9.0 

0.99567 

0.798 

76.58  b 

40 

1.00 

0.506 

1.2568 

7.039 

41.68 

0.85 

4.388 

1.2483 

5.240 

43.34 

0.75 

5.7 

1.2353 

4.743 

45.49  a 

0.50 

8.4 

1.2082 

3.212 

49.64  a 

0.25 

12.00 

1.1579 

1.928 

56.83 

0.00 

10.0 

0.99224 

0.653 

73.15  b 

a Reference  67 

b Reference  73 


The  values  for  sulfolane  at  30°C  are  in  accord  with  those 

9  10  74 
found  in  the  literature.  '   '     At  40°C,  the  dielectric 

constant  is  in  good  agreement .   The  viscosity  is  found  to 
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be  somewhat  low  and  the  density  somewhat  high  (7.953  cP 

-1  74 

and  1.2534  g  mL  "  respectively)    compared  to  the  values 

available  from  the  literature.   The  experimental  values 

were  those  used  for  the  calculations. 

It  can  be  seen  that  as  the  mole  fraction  of  water 

increases,  the  viscosity  and  density  decrease  and  the 

dielectric  constant  increases.   This  is  what  one  would 

expect.   However,  the  decrease  in  density  with  increasing 

water  content  is  not  as  great  as  would  be  expected  were  it 

simply  a  case  of  mixing  the  two  components.   There,  the 

density  would  be  expected  to  be  the  linear  combination  of 

the  two  component  densities: 

p  =  X  p   +  X  p  (53) 

S  S      WW 

where  s  and  w  denote  sulfolane  and  water  respectively.   The 
density  of  the  solvent  mixture  vs.  mole  fraction  of  sul- 
folane is  plotted  in  Figure  11a.   The  deviation  from 
ideality,  called  the  "excess  density"  is  plotted  in  Figure 
lib.   This  excess  density  is  defined  by  the  equation: 

p   =  p   -  (X  p   +Xp)  (54) 

Ke    pm    v  s*s    wKw 

where  p   is  the  measured  density, 
m 

It  will  be  noted  that  the  maximum  in  the  plot  of 
excess  density  occurs  at  about  33  mole-%  sulfolane,  at  all 
of  the  temperatures  measured.   The  deviation  from  ideality 
has  been  taken  to  indicate  association  between  the  two 
components  of  the  system.     Geddes    has  indicated  that 
the  composition  of  the  system  at  the  point  of  highest 


Figure  11a 

A  plot  of  the  density  of  the  sulf olane-water  solvents  vs. 
the  mole  fraction  of  the  sulfolane  contained  therein  at 
various  temperatures.   The  plot  for  the  densities  at  25°C 
is  not  distinguishable  from  that  for  30°C  with  the  scale 
employed. 


Figure  lib 

A  plot  of  the  excess  density  of  the  sulfolane-water  solvents 
vs.  the  mole  fraction  of  the  sulfolane  contained  therein  at 
various  temperatures.   See  text  for  definition  of  excess 
density . 
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deviation  from  ideality  is  that  of  the  associated  species. 
In  this  case,  that  would  correspond  to  sulfolane- 2H„0. 

Figure  12a  is  a  plot  of  the  dielectric  constant  vs. 
the  mole  fraction  for  the  sulfolane-water  solvent  system. 
The  excess  dielectric  constant ,  defined  in  an  analogous 
manner  to  the  excess  density,  is  plotted  in  Figure  12b. 
Here  the  deviation  is  negative  and  again  has  a  maximum 

deviation  at  about  33  mole-%  sulfolane.   This  is  in  accord 

77 
with  the  results  of  Reynaud,    and  gives  further  validity 

to  the  assumption  of  an  associated  species  with  the  com- 
position sulfolane* 2H?0. 

Figures  13a  and  13b  are  plots  of,  respectively,  the 
viscosity  and  excess  viscosity  of  the  sulfolane-water 
solvent  system.   Here  the  situation  is  more  confused  than  in 
the  previous  two  cases.   There  is  a  maximum  deviation,  but 
this  shifts  with  the  temperature.   Although  drawing  con- 
clusions from  this  plot  is  less  certain  than  with  the  plots 
of  excess  density  and  excess  dielectric  constant ,  one  thing 
stands  out.   The  addition  of  just  a  small  amount  of  water 
causes,  at  30°C,  a  drastic  decrease  in  the  viscosity  of  the 
mixture,  indicating  a  decrease  in  the  structure  of  the 
solvent.   This  is  true  to  a  lesser  extent  at  40°C.   It  has 

been  shown  that  sulfolane  is  an  excellent  "structure- 

78 
breaker"  with  regard  to  water;    apparently  the  same  can  be 

said  for  water  with  regard  to  sulfolane.   It  has  been  shown 

as  well  that  the  heat  of  mixing  for  sulfolane  and  water  is 

endothermic ;  this  would  correspond  to  heat  being  extracted 


Figure  12a 

A  plot  of  the  dielectric  constant  of  the  sulfolane-water 
solvents  vs.  the  mole  fraction  of  sulfolane  contained 
therein  at  various  temperatures.   The  25°C  plot  terminates 
at  X  =  0.85  since  sulfolane  (X  =  1.0)  is  a  solid  at  this 
temperature. 


Figure  12b 

A  plot  of  the  excess  dielectric  constant  of  the  sulfolane- 
water  solvents  vs.  the  mole  fraction  of  sulfolane  contained 
therein  at  various  temperatures.   See  text  for  definition 
of  excess  dielectric  constant. 
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Figure  13a 

A  plot  of  the  viscosity  of  the  sulf olane-water  solvents 
vs.  the  mole  fraction  of  the  sulfolane  contained  therein 
at  various  temperatures.   The  25°C  plot  terminates  at  X 
=  0.85  since  sulfolane  (X  =  1.0)  is  a  solid  at  this 
temperature . 


Figure  13b 

A  plot  of  the  excess  viscosity  of  the  sulf olane-water 

solvents  vs.  the  mole  fraction  of  the  sulfolane  contained 

therein  at  various  temperatures.  See  text  for  definition 
of  excess  viscosity. 
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from  the  solution  to  break  the  hydrogen-bonded  water 
structure.   As  will  be  shown  in  the  next  section,  this  is 
significant.   It  indicates  that  while  there  is  interaction 
between  the  sulfolane  and  water,  it  is  a  structure-breaking 
process.   For  other  dipolar-aprotic  solvents,  such  as  NM„Py 
which  is  discussed  next,  just  the  reverse  is  true. 

Pure  sulfolane  is  a  somewhat  inconvenient  solvent  to 

9 
work  with.   Since  it  melts  at  28.5°C,   most  work  has  been 

done  at  30°C,  the  first  even  temperature  at  which  sulfolane 

is  a  liquid.   However,  the  drastic  decrease  in  viscosity 

with  temperature  (drastic  when  compared  to  sulf olane-water 

solvents)  seems  to  indicate  that,  at  30°C,  the  sulfolane 

still  retains  a  great  deal  of  the  structure  that  it  had  as 

a  solid.   This  will  be  explored  further  when  the  conductance 

results  are  discussed. 

NM„Py  System 

l-Methyl-2-pyrrolidinone  (also  N-methyl-2-pyrrolidinone , 
abbreviated  as  NM„Py  herein)  is  a  dipolar  aprotic  solvent 
similar  in  many  respects  to  sulfolane.   Both  solvents  have 
five  member  heterocyclic  ring  structures;  both  have  double- 
bonded  oxygens.   The  structures  appear  in  Figure  14. 

The  physical  properties  of  NM2Py  and  NM,,Py-water 
solvents  are  summarized  in  Table  VI.   For  the  most  part, 
these  are  experimental  values,  supplemented  with  literature 
values  where  appropriate. 
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Figure  14 
Structures  of  sulfolane  and  NM0Py . 


TABLE  VI 
Physical  properties  of  NM„Py-water  solvents  at  25°C 


NM2Py 

K 

(S  cm-1) 

Density 
(g  mL"  ) 

Viscosity 

Dielectric 
Constant 

(cP) 

(  x  10"7) 

1.00 

3.637 

1.0277 

1.648 

32.06 

0.75 

7.759 

1.0348 

2.269 

37.12 

0.50 

13.86 

1.0432 

3.711 

45.01 

0.33 

28.39 

1.0477 

5.000 

52.25 

0.25 

12.12 

1.0451 

4.843 

57.00 

0.00 

7.0 

0.99707 

0.890 

78.36   a 

a Reference  73 

The  experimental  values  for  the  physical  properties  of 

pure  NM„Py  are   in  excellent  agreement  with  those  reported 

•  +u      i  -4-    *.     77,79 
in  the  literature.   ' 

It  can  be  seen  that,  as  the  mole  fraction  of  water 

present  in  the  solvent  increases,  both  the  density  and  the 

viscosity  increase  and  eventually  pass  through  a  maximum  at 

33  mole-%  NM„Py.   Figure  15a  is  a  plot  of  the  density  vs. 

the  mole  fraction  of  NMpPy  for  this  solvent  system,  while 


Figure  15a 

A  plot  of  the  density  of  the  Ni^Py-water  solvents  vs.  the 
mole  fraction  of  Nft^Py  contained  therein  at  25°C. 


Figure  15b 

A  plot  of  the  excess  density  of  the  NiV^Py-water  solvents 
vs.  the  mole  fraction  of  Nt^Py  contained  therein  at  25°C. 
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Figure  15b   is  a  plot  of  the  excess  density,  as  previously 
defined.   Figures  16a  and  16b  are,  respectively,  plots  of  the 
viscosity  and  excess  viscosity  for  these  systems.   Both 
excess  functions  pass  through  a  sharp  maximum  at  33  M% 
NMpPy ,  indicating  a  structure  for  the  associated  species  of 
NMpPy  •  2Hr,0,  analogous  to  that  found  for  the  sulfolane 
systems.   As  can  be  seen  in  Figures  17a  and  17b  (plots  of, 
respectively,  the  dielectric  constant  and  the  excess 
dielectric  constant),  this  excess  function  also  reaches  a 
maximum  value  at  33  M%.   This  finding  is  in  accord  with  the 

density  and  viscosity  results,  and  with  the  results  of 

13 

Assarson  and  Eirich,    who  demonstrated  through  heat  of 

mixing  and  phase  diagram  studies  the  formation  of  the 
aforementioned  species  NM„Py • 2H?0.   Therefore,  it  is 

believed  that  these  results  for  the  dielectric  constants 

77 
must  supersede  those  reported  by  Reynaud,    who  showed  a 

maximum  for  the  excess  dielectric  constant  at  X„,.  ~   =0.2, 

NM„Py       ' 

corresponding  to  the  species  NM2Py  •  4HpO. 

Although  the  species  formed . apparently  have  the  same 
stoichiometry ,  the  associative  processes  for  NIVUPy-HgO 
and  sulf olane-H90  must  be  quite  different.   As  was 
mentioned  in  the  previous  section,  it  is  believed  that  the 
association  of  sulfolane  with  water  is  a  structure-breaking 

process.   When  mixing  NM„Py  and  water,  it  was  quite  evident 

78 
that  the  solution  became  warm.   It  has  been  reported  '  that 

the  heat  of  solution  of  NM„Py  in  water  is  highly  exothermic, 


Figure  16a 


-  & 


A  plot  of  the  viscosity  of  the  NiV^Py-water  solvents  vs. 
the  mole  fraction  of  NJ^Py  contained  therein  at  25°C. 


Figure  16b 


A  plot  of  the  excess  viscosity  of  the  NlV^Py-water  solvents 
vs.  the  mole  fraction  of  NA^Py  contained  therein  at  25°C. 
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Figure  17a 

A  plot  of  the  dielectric  constant  of  the  Ni^Py-water 
solvents  vs.  the  mole  fraction  of  NN^Py  contained  therein 
at  25°C. 


Figure  17b 

A  plot  of  the  excess  dielectric  constant  of  the  Ni^Py-water 
solvents  vs.  the  mole  fraction  of  Ni^Py  contained  therein 
at  25°C. 
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consistent  with  the  empirical  observation.   This,  coupled 
with  the  large  increase  in  solution  viscosity,  makes  it 
evident  that  the  association  between  Nlv^Py  and  water  is 
a  structure-making  process,  probably  involving  hydrogen 
bonding  at  the  carbonyl  oxygen. 

Pure  NM2Py  has  been  shown  to  exist  primarily  as  a 

SO 
monomer.     As  is  evident  from  the  viscosity  curve,  addition 

of  small  amounts  of  water  causes  a  tremendous  increase  in 

the  structure  of  the  solvent,  indicating  a  greater  tendency 

for  the  NM„Py  to  associate  with  water  rather  than  itself. 

The  effect  of  the  various  solvent  properties  which  have 

been  discussed,  such  as  association  between  the  two 

components  of  the  mixture,  is  discussed  in  a  following 

section  of  this  chapter. 


-87- 
Conductance  Curves 

When  conductance  data  are  presented  graphically  it 
is  customary  to  plot  the  equivalent  conductance,  A,  vs. 
the  square  root  of  the  concentration.   For  a  dilute  solution 
of  a  completely  dissociated  electrolyte,  this  would  have 
a  form  very  close  to  the  Debye-Huckel-Onsager  equation, 


o 


A  =  Aw  -  Sc2  (21) 

with  the  slope  of  the  line  being  -S  and  the  intercept  being 
the  limiting  equivalent  conductance,  A  .   Experimentally, 
electrolytes  that  have  negligible  association  give  straight- 
line  conductance  plots  which  extrapolate  quite  closely  to 
the  actual  A   value.   Such  linear  plots  generally  lie 
above  the  Debye-Hiickel-Onsager  slope. 

In  their  book,  "Electrolytic  Conductance,"  Fuoss  and 

46 
Accascina   have  described  the  behavior  of  the  conductance 

curve  as  the  association  constant,  K  ,  is  increased.   As  K 

'   a  a 

increases  from  the  negligible  case  described  above,  the 
conductance  curve  falls  below  the  D-H-0  slope   and  acquires 
an  increasingly  negative  deviation  as  K   increases.   They 
have  divided  the  curves  into  six  types,  making  very  careful 
mention  of  the  fact  that  there  is  a  continuous  transition 
between  the  various  types  of  curves.   Data  from  this  study 
to  fit  each  type  but  one  of  the  curves  described  by  them 
were  available  and  they  are  presented  here. 

Figure  18  is  an  example  of  the  Type  I  conductance  curve 
It  is  for  KC1  in  water  at  25°C,  and  it  is  representative  of 
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the  curves  for  strongly  dissociated  electrolytes.   At  very 
low  concentrations,  the  curve  lies  on  the  D-H-0  slope.   As 
the  concentration  is  increased,  a  positive  deviation  is 
noted,  indicative  of  the  dissociated  electrolyte. 

Figure  19  represents  the  Type  II  curve.   Here,  K   is 
no  longer  negligible  (characteristic  of  solvents  with 
median  values  for  the  dielectric  constant)  and  the  curve 
lies  slightly  below  the  D-H-0  slope.   Extrapolation  to 
zero  concentration  still  yields  a  good  value  for  A  , 
however.   The  curve  presented  is  for  HC1  in  25  M%  sulfolane- 
water,  at  30°C,  and  has  an  experimentally  determined  K   of 
2.37  mol    L.   The  curve  is  typical  for  electrolytes  with 
association  constants  in  this  range,  where  ion  association 
is  just  starting  to  take  place. 

Figure  20  represents  a  Type  III  curve,  in  this  case 

for  HC1  in  50  M%  sulf olane-water  at  30°C.   The  experimentally 

determined  K   for  this  curve  is  27.73  mol    L,  and  this 
a 

Type  III  behavior  will  generally  be  seen  with  association 
constants  ranging  between  15  and  200  mol    L.   Here 
extrapolation  of  the  experimental  curve  will  yield  a  value 
for  A   which  is  high  compared  to  the  actual  value.   It  is 
at  this  point,  where  K   has  become  sufficiently  large  to 
influence  the  extrapolation,  that  a,  the  degree  of 
dissociation,  should  be  incorporated  into  the  theoretical 
treatment  used  to  describe  the  behavior  of  the  curve. 

Figure  21  is  an  intermediate  curve,  representing  the 
transition  between  Types  III  and  IV.   This  curve  is  for  HC1 
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in  75  M%  sulfolane-water  at  30°C,  with  K   determined  to 

a 

be  288.13  mol    L.   In  this  type  of  behavior,  the  inflection 

point  in  the  curve  has  started  to  move  into  the  low 

concentration  region,  and  extrapolation  of  the  curve  to 

zero  concentration  will  yield  a  value  for  A   that  is  rather 

higher  than  the  actual  value. 

Figure  22  is  Type  IV  behavior,  characterized  by  the 

inflection  point  of  the  curve  being  in  the  low  concentration 

region.   Extrapolation  to  zero  concentration  gives  a  value 

for  the  limiting  conductance  that  is  substantially  in  error. 

This  curve  is  for  HC1  in  85  M%  sulfolane-water  at  40°C, 

with  a  K   of  1336.2  mol-1  L. 
a 

For  all  of  the  curves  described  above,  the  full  Pitts 
equation  was  found  to  yield  excellent  agreement  between  the 
calculated  and  the  experimental  values  for  the  equivalent 
conductance.   The  expanded  Pitts  equation  and  the  expanded 
Fuoss-Hsia  equation  were  also  found  to  yield  satisfactory 
agreement.   With  the  computer  equipment  available  today, 
utilizing  the  proper  program  allows  the  use  of  a  single 
theoretical  treatment  for  conductance  behavior,  rather  than 
the  several  treatments  recommended  by  Fuoss  and  Accascina 
in  their  treatment  of  the  conductance  curves. 

Figure  23  is  representative  of  Type  V  behavior.   This 
curve  is  for  HC1  in  sulfolane  at  40°C,  and  is  characteristic 
of  Type  V  behavior,  showing  very  small  values  for  A,  and 
very  high  curvature  in  the  low  concentration  region.   Fuoss 
and  Accascina  state  that  there  should  also  be  a  minimum 
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present ,  and  attribute  this  behavior  to  a  very  high 
association  constant  coupled  with  triple  ion  formation. 
That  is  not  noted  here,  however,  and  as  will  be  discussed 
in  a  following  section,  triple  ion  formation  in  this  case 
is  unlikely. 

Figure  24  is  another  example  of  Type  V  behavior,  in 
this  case  for  HC1  in  sulfolane  at  30°C.   Here,  the  minimum 
due  to  triple  ion  formation  is  noticeable;  as  the  concen- 
tration increases,  the  concentration  of  triple  ions 
increases,  thus  increasing  the  conductance.   Neither  the 
program  for  the  full  Pitts  equation  nor  that  for  the 
expanded  equations  was  capable  of  predicting  this  Type 

V  behavior.   Rather,  a  differing  extrapolative  method  was 

42 
utilized,  that  of  Shedlovsky  and  Kay,    and  will  be 

discussed  in  a  following  section. 

Type  VI  is  associated  with  solvents  of  the  hydrocarbon 

type,  with  very  low  dielectric  constants,  and  is 

characterized  by  increasing  conductance  with  concentration 

coupled  with  a  complex  curve  form. 
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Ion  Size  Parameter 

General 

In  the  preceding  chapters,  mention  has  been  made  of 
the  so-called  ion  size  parameter.   This  is,  in  fact,  a 
parameter  introduced  into  the  extended  Debye-Hiickel 
expression  for  the  activity  coefficient  and  into  the  main 
conductance  equation  via  kcl   to  account  for  the  distance  of 
closest  approach  of  the  ions  in  solution.   There  is  a  large 
element  of  arbitrariness  in  the  choice  of  a  value  for  this 
parameter,  and  in  the  course  of  this  study  several  methods 
were  employed  to  obtain  a  reasonable  value. 

A  change  in  the  value  of  a.   will  cause  a  minor  variation 
in  the  value  of  A  ,  and  a  greater  change  in  the  value  of 
K  ,  produced  by  application  of  the  conductance  equations 
to  the  data.   For  electrolytes  in  high  dielectric  constant 
media,  such  as  water,  it  has  been  noticed  that  as  the  ion 
size  parameter  is  varied,  the  standard  deviation  between 
the  calculated  and  experimental  conductances,  a,,  goes 
through  a  minimum.   This  minimum  in  a,    is  taken  as  a 
criterion  of  best  fit ,   '    and  the  value  of  the  ion  size 
parameter  at  that  point  is  used  to  generate  the  calculated 
conductance  parameters. 

This  procedure  was  used  in  this  study  for  Tris-HCl 
in  water.   A  value  of  a   was  selected,  and  the  conductance 
parameters  were  generated.   The  value  of  a   was  incremented, 
and  the  procedure  was  repeated.   The  value  of  a  which 
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yielded  a  minimum  a.  was  the  one  selected  for  use  in 
generating  the  results. 

When  dealing  with  solvents  possessing  medium  to  low 
dielectric  constants,  however,  things  become  more  complex. 
It  has  been  noted  that  varying  a   will  yield   a  minimum  in 

o,,    but  it  will  be  extremely  shallow,  or  even  a  straight 

82 
line.     In  other  cases,  no  minimum  is  observed;  rather, 

there  is  a  maximum!   In  such  cases,  the  choice  of  a  is 

rather  arbitrary,  with  the  value  selected  being  "reasonable," 

whatever  that  means. 

An  example  of  this  behavior  of  a.  with  respect  to  the 

value  of  a   follows.   For  HC1  in  75  M%  sulfolane-water  at 

—  8 
25°C,  as  shown  in  Table  I,  with  an  ion  size  of  5.99  x  10 

cm,  a.  =  0.121.   If  this  value  for  the  ion  size  is  changed 

—8 
to  3  x  10   cm,  a.  =  0.125.   Very  clearly,  varying  the 

value  of  the  ion  size  parameter  is  not  a  suitable  method 

for  obtaining  a  reasonable  value,  in  this  case. 

Because  of  the  problem  in  choosing  a  reasonable  value 

for  the  ion  size  parameter,  several  authors     '  have 

proposed  using  the  Bjerrum  critical  distance  as  the  ion 

size  parameter.   That  is  the  procedure  used,  for  the  most 

part,  in  this  study. 

Bjerrum' s  distance 

A  fundamental  question  when  dealing  with  conductance 
and  ion  association  which  must  be  answered  is,  what  is  an 
ion  pair?   At  what  point  are  ions  considered  paired,  since 
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at  that  point  they  will  cease  to  contribute  to  conductance? 
Must  they  be  in  contact,  or  merely  nearby?   Most  of  the 

work  which  has  been  done  with  respect  to  a  solution  to  this 

84 
problem  has  been  based  on  the  early  work  of  Bjerrum. 

It  was  Bjerrum  who  proposed  that  all  oppositely  charged 

ions  within  a  certain  critical  distance  be  considered 

paired  and  no  longer  contributing  to  conductance.   Basing  his 

work  upon  a  Boltzmann  distribution  function  for 

electrostatic  forces,  Bjerrum  obtained  a  probability 

function  for  finding  an  ion  at  a  distance  q  from  a  central 

ion.   A  minimum  was  noted  in  this  function,  corresponding, 

for  a  1:1  electrolyte,  to: 

2  2 
q  -  f^fr  (55) 

Ions  closer  together  than  this  critical  distance  are,  in  the 
Bjerrum  treatment,  considered  to  be  associated,  since  the 
probability  of  association  increases  very  rapidly  below  this 
distance. 

Applying  the  law  of  mass  action  to  free  and  associated 
ions,  Bjerrum  derived  an  equation  for  the  degree  of 
association,  1  -  a: 

1  "  a  =  1000  J  exp("  IkT7)r  dr  (°6) 

a 

From  this,  an  expression  can  be  derived  for  the  association 

constant.   This  expression  will  be  based  totally  on 

electrostatic  considerations: 
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4-rrNa   e 
a    1000   b  ^°n 

where : 

b  =H  (58) 

85 
Fuoss,   however,  has  stated  that  only  those  ions  which 

are  actually  in  contact  with  one  another  can  be  considered 

to  be  associated.   His  analysis  of  the  problem  yields  the 

equation : 

4ttNo.  e 
K  =  ?  (59') 

a     3000  ^Dy; 

It  has  been  pointed  out86  that  these  electrostatic 
models,  in  all  but  the  simplest  cases,  are  totally  in- 
adequate for  representing  real-life  behavior.   Indeed,  this 
was  found  to  be  the  case  in  this  study  and  will  be  further 
discussed  when  the  results  are  examined.   In  any  case, 

quoted  results  for  A   and  K   should  always  make  mention 

a 

of  the  ion  size  value  used  in  the  calculations,  and  that  rule 
is  adhered  to  in  this  study. 
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Choice  of  Conductance  Equation 

After  all  the  data  are  gathered,  and  after  a  method 
of  arriving  at  a  reasonable  value  for  the  ion  size  parameter 
is  decided  upon,  there  still  remains  the  choice  of  which 
conductance  equation  to  use.   Each  has  advantages  and  disad- 
vantages.  The  expanded  equations,  both  Pitts  and  Fuoss-Hsia, 
take  very  little  time  on  the  computer.   However,  they  do 
tend  to  break  down  at  high  values  of  kcl.      The  full  Pitts 
equation  is  generally  applicable  to  conductance  data  but 
takes  a  rather  long  time  before  the  results  are  computed. 
Of  course,  where  the  expanded  equations  work,  they  could  be 
employed,  and  where  they  do  not,  the  full  Pitts  equation 
program  could  be  utilized. 

Perhaps  the  best  method  of  deciding  which  equation 
to  use  is  to  examine  how  well  the  equation  can  fit  the 
data;  that  is,  which  equation  can  theoretically  synthesize 
equivalent  conductance  values  closest  to  those  experimen- 
tally obtained.   In  this  way,  the  logical  choice  of 
equation  would  be  that  which  gives  the  best  fit. 

Table  VII  presents  the  results  for  HC1  in  sulfolane- 
water  solvents,  calculated  by  all  three  equations. 
Examination  of  the  values  of  a.  reveals  that,  in  general, 
the  full  Pitts  equation  gives  a  more  satisfactory  fit  to 
the  conductance  data  than  either  of  the  expanded  equations. 
This,  coupled  to  the  fact  that  the  expanded  equations  broke 
down  for  several  of  the  solvent  systems  studied,  is  the 
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rationale  by  which  the  full  Pitts  equation  was  chosen  for 
the  analysis  of  data  in  this  study.   All  results  presented, 
with  the  exception  of  those  labeled  otherwise,  such  as  in 
Table  VII,  are  generated  with  the  full  Pitts  equation, 
using  the  Bjerrum  critical  distance  as  the  ion  size 
parameter. 
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2-Methoxyethanol  -  Water  Solvent  System 


General 

The  results  for  this  study  were  presented  in  Table 
III.   Those  results  utilize  the  Bjerrum  critical  distance, 

which  in  the  solvent  studied,  80  W%  MetOH-HgO,  is 

—  8 
8.7  x  10  '  cm.   Analysis  of  the  data  by  varying  the  ion 

size  parameter  showed  a  minimum  in  the  plot  of  a.    vs.  a. 

As  an  example  of  how  the  varying  value  of  the  ion  size 

parameter  affects  the  results ,  the  conductance  parameters 

generated  by  this  method  are  listed  below. 

TABLE  VIII 

Calculated  conductance  parameters  for  the  80  W%  MetOR-HoO 
solvent  system,  utilizing  the  variation  in  a   value  method. 

Solute  A°  K  a  a. 

(S  cm   mol   )     (mol    L)    (cm  x  10   ) 

HC1  53.066 

NaCl         31.651 

Sodium  acetate  was  found  to  give  a  maximum  in  a.  rather 

than  a  minimum. 

Variation  of  the  ion  size  parameter  from  the  Bjerrum 

—  8 
distance  of  8.7  x  10  '  cm  to  about  half  that  value  has 

caused  a  change  of  only  about  0.5%  in  A  ;  at  the  same  time, 

the  values  for  K   have  changed  by  almost  100%.   The  Bjerrum 

distance  is  used  throughout  this  study,  and  the  results  are 

those  found  in  Table  III. 
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10.28 

4.5 
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The  value  for  the  association  constant  of  HC1 , 
14.81  mol    L,  corresponds  to  a  pK,  of  1.17.   This  is  in 

"17  1  8 

good  agreement  with  the  results  of  Merken  et  al. ,   '   '  who 

report  a  pK ,  of  1.20.   Analysis  of  their  conductance  data 

with  the  full  Pitts  equation  yielded  value  for  pK,  of  1.22, 

still  in  excellent  agreement. 

It  is  of  interest  to  note  that,  in  a  previous  paper, 

Thun  and  his  coworkers    estimated  K   for  HC1  from  EMF 

a 

measurements.   Their  value  of  4.902  mol-   L  is  somewhat 

lower  than  that  found  conductimetrically ,  pointing  out 

the  value  of  direct  measurement  over  roundabout  estimates. 

It  is  possible,  using  the  Bjerrum  and  Fuoss  theories 

of  association,  to  obtain  a  value  for  K  based  solely  on 

a 

electrostatics.   For  the  80  W%  MetOH-H?0  system,  using  the 

Bjerrum  distance,  this  is  18.40  mol    L  with  the  Bjerrum 

theory  and  12.27  mol    L  with  the  Fuoss  theory.   It  can  be 

seen  that  there  is  fair  agreement  with  the  experimental 

result,  indicating  that  electrostatic  effects  go  a  long  way 

toward  accounting  for  the  association  of  HC1  in  the  system. 

It  has  been  postulated    that  in  this  system  the  proton 

is  highly  solvated  by  water.   If  that  is  the  case,  the  proton 

and  the  CI   ion  would  be  held  further  apart,  and  may  explain 

why  the  association  constant  for  HC1  is  less  than  for  NaCl . 

The  value  for  K  of  HAc  of  1.24  x  108  mol~   L  was  not 

a 

derived  by  the  full  Pitts  equation.   This  electrolyte  showed 
Type  V  behavior,  as  explained  in  a  preceding  section,  and 
the  Pitts  equation  was  unable  to  account  for  this  behavior. 
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Use  was  made  here  of  the  extrapolative  method  devised  by 

42 

Shedlovsky  and  Kay    for  the  treatment  of  their  data  for 

HAc  in  methanol-water  solvents,  and  it  is  applicable  to 
systems  where  the  conductance  due  to  the  solvent  is  an 
appreciable  fraction  of  the  total  conductance. 


Shedlovsky-Kay  Extrapolation 

This  extrapolative  method  for  the  determination  of 

K   for  very  weak  acids  is  quite  useful,  provided  one  major 
a 

condition  can  be  fulfilled.   For  application  of  this  method 
it  is  necessary  to  have  an  estimate  of  A  .   The  usual 
method  of  obtaining  this  is  to  use  the  Kohlrausch  law  of 
independent  migration.   According  to  this,  the  moving  ions 
do  not  influence  each  other  and,  consequently,  the  limiting 
conductance  of  the  ions  remains  constant  despite  the  other 
ions  in  solution.  Thus  A   can  be  determined  by  combining 


me 


asurements  on  other  electrolytes.   For  80  W%  MetOH-H„0, 


the  value  of  A   for  HAc  was  synthesized  from  the  values 
for  HC1,  NaCl ,  and  NaAc ,  as  follows. 

The  limiting  conductance  of  HAc  is  made  up  of  the 
contributions  of  the  ionic  conductances: 

A^a   =  ^   +  *a  (60) 

HAc     H     Ac 

We  can  obtain  this  combination  of  ionic  conductances  by 

the  combination  of  A   values  for  the  aforementioned 

electrolytes : 

AHAc  '  AHC1    ANaAc  "  ANaCl  "  AH  +  AC1    ANa    AAc  " 

(61) 
ANa  "  AC1  _  AH    AAc 
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For  this  solvent  system,  the  procedure  yields  a  value  of 

O  "1 

A   for  HAc  of  47.064  S  cm   mol   .   Using  this  value,  the 
Shedlovsky-Kay  extrapolation  can  be  applied. 

In  this  extrapolation,  the  equivalent  conductance 
which  contains  the  solvent  correction  is  not  used,  but 
rather  one  thousand  times  the  specific  conductance.   This 
quantity,  abbreviated  L,  is  plotted  vs.  a  function  of  the 
square  root  of  the  concentration.   It  is  during  this 
extrapolation,  in  calculating  the  value  of  the  function  of 
c2 ,  that  the  solvent  correction  is  applied.   With  L   as 
one  thousand  times  the  specific  conductance  of  the  solvent, 
the  Shedlovsky-Kay  extrapolation  equation  becomes: 

A°K  *c*       *Q      L    1 
L-— £^(1-^1-  ■£))*-  LQ  (62) 

where  f  is  the  activity  coefficient,  A   is  L  divided  by 
the  concentration,  and  S  is  a  function  of  c  and  A  . 

The  computer  program,  "SHEDKAY,"  was  designed  to 
follow  the  recommendations  of  the  authors  of  the  paper  as 

to  the  proper  extrapolation  technique.   Using  that  program, 

8     —1 
a  value  was  obtained  for  K   of  1.24  x  10   mol    L,  or  a 

a 

pK,  of  8.09.   This  may  be  compared  to  values  of  6.71  in  80  W% 
methanol-water    and  6.98  in  80  W%  ethanol-water .        It 
can  be  seen  that  in  all  three  alcohols  acetic  acid  exists 
in  an  essentially  ion-paired  form,  and  in  all  three  it  is 
far  more  associated  than  in  water  (pK.,  =  4.76).   This  is  in 
keeping  with  the  idea  that  water  can  solvate  the  proton  more 
easily  than  can  the  bulky  alcohols.   The  A   value  in  SO  W% 
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MetOH-H^O  is  much  higher  than  would  be  expected  with  an 

essentially  associated  eletrolyte,  and  may  be  attributed 

17  IF 
to  the  limited  proton  jump  mechanism  in  the  alcohols.   ' 

A  preliminary  value  for  pK,  of  HAc  in  80  W%  Met0H-Ho0, 

71 
determined  via  EMF  studies  in  this  laboratory    is  7.39,  a 

considerable  difference  from  that  determined  by  this  study. 

The  reason  for  the  difference  remains  unexplained. 

At  25°C,  the  K   for  Tris-HCl  in  water  was  found  to  be 
'       a 

1.01  mol-1  L,  and  in  80  W%  MetOH-H90  to  be  29.46  mol-1  L. 
These  results  may  reflect  the  lesser  ability  of  the  alcohol 
to  solvate  the  TrisH   ion  which,  with  three  OH  groups,  ca.n 
enter  into  hydrogen  bonding  with  the  water  structure.   This 
interaction  with  the  water  structure  is  seen  in  the  low 

limiting  ionic  conductance  found  for  the  TrisH   ion  in 

72 

water.     It  is  interesting  to  note  that,  in  water  at  25°C, 

using  the  ion  size  determined  for  Tris-HCl,  the  Bjerrum 

association  constant  is  0.78  mol    L  and  the  Fuoss 

association  constant  is  0.70  mol-   L.   Both  are  reasonably 

close  to  the  experimentally  determined  K   of  1.01  mol    L, 

a 

indicating  that  the  association  is  primarily  electrostatic 

in  nature.   This  is  not  the  case  in  80  W%  MetOH-H90,  where 

Bjerrum  and  Fuoss  K  values  are,  respectively,  18.4  and 

a 

12.3  mol-   L.   The  experimentally  determined  constant  of 
29.46  mol    L  indicates  that  electrostatics  cannot  wholly 
explain  the  association  found.   More  will  be  said  about  this 
later o 


■117- 


Sulfolane-Water  Solvent  System 

General 

Sulfolane  is  probably  the  perfect  example  of  a  dipolar 
aprotic  solvent.   A  number  of  conductimetric  studies  have 

been  done  using  sulfolane  as  the  solvent ,  almost  universally 

9  89-99 
at  30°C.  '         Very  little  work  has  been  done  with 

-,-p  i       +.     -^-14.   95-98,100,101   T17.^u  ... 
sulf olane-water  mixed  solvents.  '      With  this  in 

mind,  a  study  of  the  conductance  of  HC1  in  sulfolane  and 
sulf olane-water  solvents  at  several  temperatures  was  under- 
taken . 

The  conductance  results  for  this  aspect  of  the  study 
are  summarized  in  Table  I.  Those  calculations  were  made  with 
the  Bjerrum  critical  distance  as  the  value  for  the  ion  size 
parameter  in  both  the  Pitts  equation  and  the  expression  for 
the  activity  coefficient.   The  results  for  pure  water  were 
obtained  by  applying  the  full  Pitts  equation  to  the  data  of 
Owen  and  Sweeton  '  and  Stokes.   '     HC1  in  pure  sulfolane 
exhibited  Type  V  conductimetric  behavior  and  is  discussed 
separately . 

Figure  25  shows  the  behavior  of  the  limiting  equivalent 
conductance  with  respect  to  the  mole  fraction  of  sulfolane 
in  the  solvent.   As  can  be  seen,  the  addition  of  small 
amounts  of  sulfolane  to  water  cause  a  dramatic  decrease  in 
the  limiting  conductance  of  HC1.   This  may  be  attributed  to 
the  disruption,  by  sulfolane,  of  the  water  structure,  and  a 
consequent  dramatic  decrease  in  the  probability  of  the  proton 
jump  mechanism  taking  place. 
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It  is  the  proton  jump  mechanism  that  is  responsible 
for  the  anomalously  high  conductances  of  the  proton  and 
the  hydroxide  ion  in  water.   This  is  due  to  the  fact  that 
these  ions  are  actually  fragments  of  the  solvent,  water. 
With  the  hydrogen-bonded  structure  of  water,  there  is  no 
need  for  the  ions  to  move.   Rather,  the  proton  (or 
hydroxide  ion)  is  transferred  from  a  hydronium  ion  to  a 
water  molecule,  which  transfers  it  on,  etc.   The  result  is 
no  great  movement  of  the  ion  per  se ,  but  a  very  rapid 
movement  of  the  charge  along  the  chain  of  hydrogen-bonded 
water  molecules.   The  proton  jump  mechanism  is  diagrammed 
in  Figure  26. 

H    H        H    H 
H0H0H0H0H  t   H0H0H0H0H 
H    H        H    H 

Figure  26 

The  proton  jump  mechanism. 

If  the  water  structure  were  disrupted,  the  length 
of  the  structural  chains  would  decrease,  and  a  decrease 
in  the  conductance  of  the  proton  would  be  expected.   This 
is  exactly  what  is  noted  when  sulfolane  is  added  to  water, 
and  ties  in  well  with  the  structural  considerations  noted 
earlier.   By  the  time  the  solvent  contains  25  M%  sulfolane, 
the  limiting  conductance  has  decreased  by  75%,    indicating 
that  the  proton  jump  mechanism  (and  hence  the  water  struc- 
ture) has  been,  for  the  most  part,  disrupted. 
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Figure  27  is  a  plot  of  the  pK  ,  (=  log  K  )  of  HC1  vs. 
the  mole  fraction  of  sulfolane.   As  can  be  seen,  there  is 
almost  a  linear  dependence  of  pK,  on  mole  fraction  of 
sulfolane  for  most  of  the  solvent  mixtures,  the  exceptions 
being  at  either  end  of  the  plot.   From  approximately  10  M% 
to  90  M%,  there  is  a  linear  dependence  on  mole  fraction  of 
sulfolane.   Since  the  amount  of  ionization  is  directly 
proportional  to  the  strength  of  solvation  of  the  ions 
involved,  it  is  quite  apparent  that  as  the  amount  of  sul- 
folane in  the  solvent  increases,  the  solvation  of  the  ions, 
in  particular  the  proton,  decreases. 

The  rather  precipitous  drop  in  pK_,  as  small  amounts 
of  water  are  added  to  sulfolane  probably  corresponds  to 

the  hydration  of  the  proton  to  form  hydronium  ions,  thus 

102 
promoting  dissociation.  Milanova  et  al      have  pointed  out 

that  as  water  is  added  to  acid  solutions  in  sulfolane,  the 

+ 


species  HCH^O)   are  formed,  with  the  primary  species  being 
H(Hr,0)4.   Once  this  species  is  formed,  the  increase  in 
dissociation  is  far  more  moderate,  probably  indicative  of 
longer  chains  of  water  molecules  forming  as  the  water  content 

is  increased. 

103 

Bates  and  Pawlak     have  pointed  out  the  linear 

dependence  of  pK,  on  the  mole  fraction  of  sulfolane  in  the 
range  X   =  0  to  0.8  for  a  variety  of  weak  acids,  similar  to 
the  dependence  found  here  for  HC1.   As  they  assumed,  the 
pp[   does  rise  significantly  (to  say  the  least)  between 
mole  fractions  of  0.8  and  1.0. 
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96  97 
Petrella  and  his  co-workers   '    have  examined  the 

conductance  of  HC1  in  sulf olane-water  mixtures  at  35°C, 
which  would  complement  the  25,  30,  and  40°C  measurements 
of  this  study.   Unfortunately,  they  avoided  systems  where 
association  was  not  negligible  and,  as  a  result,  inter- 
comparison  of  results  is  not  feasible.   Their  limiting 
conductance  values,  plotted  in  the  manner  of  Figure  25, 
fall  roughly  between  the  curves  for  30  and  40°C,  as 
expected.   However,  their  values  for  K   appear  to  be  rather 
high  compared  to  the  results  obtained  in  this  study.   For 

XQ  =  0.3116,  they  report  a  K   of  51  mol    L,  whereas  inter- 
im a 

polating  with  Figure  27  yields  a  value  of  approximately 
7.25  mol-   L.   Again,  at  X„  =  0.1583,  they  report  K   =  26 
mol    L,  and  using  Figure  27  a  value  of  less  than  1 
mol    L  is  obtained.   The  reasons  for  this  difference  may 
lie  in  the  manner  of  treating  the  raw  conductance  data, 
since  they  did  not  use  a  solvent  correction  when  cal- 
culating the  equivalent  conductances,  and  they  treated 
their  data  with   the  Fuoss-Onsager  equation  for  obtaining 
results . 

The  variation  of  K  with  temperature  for  the  sulfolane- 
rich  solvents  is  worthy  of  mention.   The  value  of  K   goes 
down  with  increasing  temperature  from  25  to  30°C,  but  then 
increases  from  30  to  40°C.   A  possible  rationale  for  this 
behavior  follows.   As  the  temperature  increases,  the 
dielectric  constant  of  the  solvent  decreases,  and  this 
will  promote  the  formation  of  ion  pairs.   At  the  same  time, 
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increasing  temperature  increases  the  Brownian  motion  within 

the  system,  which  would  tend  to  break  up  ion  pairs.   The 

two  competing  effects  may  be  responsible  for  this  strange 

behavior. 

For  these  systems,  the  electrostatic  association  models 

clearly  show  inadequacies.   At  25  M%  sulf olane-water ,  at 

25°C,  the  K   from  the  Bjerrum  theory  is  2.88  mol~   L;  from 
a 

the  Fuoss  theory  it  is  1.92  mol-   L.   The  experimentally 
determined  value  of  2.54  mol    L  is  in  rather  close 
agreement.   Interestingly  enough,  the  method  of  variation 
of  the  ion  size  parameter  mentioned  earlier  worked  well  for 
this  system.   Clearly,  the  system  may,  for  the  most  part, 
have  its  behavior  explained  in  terms  of  electrostatics. 

For  the  rest  of  the  sulf olane-water  solvents,  however, 
this  is  not  the  case.   For  50-,  75-,  and  85-M%  sulf olane- 
water  solvents  at  25°C,  the  Bjerrum  theory  predicts  values 
for  the  association  constants  of  4.30,  6.01,  and  6.63  mol  '  L 
respectively;  for  the  Fuoss  theory,  2.87,  4.01,  and  4.42 
mol    L,  respectively.   Comparison  with  the  actual 
experimental  values  shows,  quite  clearly,  that  electro- 
statics can  no  longer  explain  the  ion  association.   It  now 
becomes  more  a  matter  of  solute-solvent  interactions,  as 
previously  pointed  out. 

The  major  anomaly  on  the  plot  of  pK,  vs.  mole  fraction 
sulfolane  (Figure  27)  is  quite  clearly  at  1.0  mole  fraction 
sulfolane.   Such  a  drastic  change  in  pK,  is  evidence  of  a 
fundamental  change  within  the  solution,  and  that  is  explored 
next . 
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Pure  Sulfolane 

Since  the  first  conductance  study  in  sulfolane,  by 

89 
Burwell  and  Langford  in  1959,    there  have  been  a  number  of 

studies  of  the  solvent,  but  many  questions  yet  remain  about 

solution  processes  in  this  most  interesting  solvent. 

Most  studies  in  the  past  have  been  done  at  30°C,  since 

this  is  the  first  even  temperature  above  the  melting  point 

of  sulfolane.   Such  studies  have  shown  that  the  ability  of 

sulfolane  to  solvate  ions  is  in  direct  proportion  to  the 

size  of  the  ion.   Electrolytes  containing  large  ions,  such 

as  tetraphenylarsonium,  appear  to  be  largely  dissociated  in 

,,  ,       n  ..     9,89,91,92,98,99   .  ,.   ..    ..  .  .. 

sulfolane  solutions,  '         '        indicating  that  they 

are  solvated  well  by  the  solvent.   By  contrast,  electrolytes 

with  small  ions,  such  as  the  alkali  halides,  have  been  found 

9 
to  exist  primarily  in  undissociated  forms. 

There  have  been  attempts  to  set  up  acidity  and 

u   •  •+      i    -     -i-P-,        -+u    i  •  •+  ^  104-109 

basicity  scales  in  sulfolane,  with  limited  success. 

Examining  freezing-point  depressions  caused  by  sulfuric  acid, 

105 
Arnett     came  to  the  conclusion  that  sulfolane  will  not 

108 
protonate.   Whiting  and  his  coworkers,     in  the  course  of 

setting  up  such  a  scale,  came  to  the  conclusion  that  acids 

with  large  anions  (e.g.,  C10~   PF~ )  were  dissociating  simply, 

while  the  mineral  acids  were  subject  to  homoconjugation  as 

well : 

HC1  X   H+  +  CI"  (63a) 

HC1  +  CI"  %   HC1~  (63b) 
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Coetzee        has  used  sulfolane  as  a  medium  for  conduct- 

imetric  titrations  and  obtained  a  pK,  for  HC1  by  titration 

with  silver  perchlorate  potentiometrically .   He  obtained 

96  97 
a  value  of  14.5,   Petrella  and  coworkers   '    determined 

K   for  HC1  conductimetrically  at  35°C  and  obtained  a  value 

Si 

of  pK  ,  of  9.0.   There  is,  evidently,  some  confusion  about  the 
properties  of  sulfolane  as  a  solvent. 

The  conductance  of  sulfolane,  at  30°C,  varies  very 

2 
little  with  concentration,  remaining  at  about  0.22  S  cm 

mol   .   This  is  a  perfect  example  of  Type  V  conductance,  and 

none  of  the  theoretical  equations  was  able  to  account  for 

this  type  of  behavior.   A  different  approach  was  tried, 

that  suggested  by  Caruso  and  his  coworkers. 

The  method  of  Caruso  involves  plotting  -log  A  vs. 

-log  c.   The  intercept  of  the  plot  provides  information  about 

K  ,  and  the  slope  of  the  line  varies  with  the  nature  of  the 
a'  L 

dissociation  process.   Figure  28  shows  such  Caruso  plots 
for  HC1  in  sulfolane  at  30°C  and  40°C.   The  plot  of  40°C 
has  a  slope  of  approximately  -0.5,  indicating  a  simple 
dissociation  process  (Equation  63a) .   The  intercept  for  such 
an  ionization  process  has  the  magnitude  -§logK,A   .   The 
limiting  conductance  cannot  be  obtained  by  extrapolation 

and  is  obtained  from  literature  values.   Delia  Monica  and 

90  o  o 

coworkers    report  the  value  of  Xp,  as  9.30,  and  Ap-iQ   as 

6.68  S  cm  mol"  .   Coetzee     gives  the  limiting  conduc- 
tance  of  HC10.  as  15.   Combining  this  information  gives 

X°  =8.3  and  A°  for  HC1  =  17.6  S  cm2  mol-1.   This  yields  a 
rl 

-6     -1 
value  of  K   of  6.14  x  10    mol    L,  or  a  pK,  of  6.79. 
a  '         d 
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Since  the  process  is  a  simple  dissociation,  the 

Shedlovsky-Kay  extrapolation  for  weak  acids  was  applied. 

6     -1 
This  gave  a  value  of  2.45  x  10   mol    L  for  K  ,  or  pK,  = 

6.39.   Since  this  method  of  arriving  at  K   involves  fewer 

°  a 

approximations,  these  are  regarded  as  the  best  values  in 
this  study. 

Examination  of  Figure  28  shows  that,  at  30°C,  the  plot 
of  -log  A  vs.  -log  c  has  a  slope  of  zero.   Caruso  has 
indicated  that  this  is  typical  of  a  simple  dissociation 
coupled  with  homocon jugation  (Equations  63  a  &  b) .   This 
is  in  accord  with  the  results  of  Whiting,  mentioned  earlier. 

Benoit  and  his  coworkers  have  studied  the  system 

111—112 
HC1  in  sulfolane  at  30°C.         Utilizing  calorimetry 

and  proton  NMR ,  they  concluded  that,  in  sulfolane,  HC19 

forms,  and  that  the  formation  constant  is  approximately  350. 

Their  finding  is  in  accord  with  the  Caruso  plot  at  30°C. 

For  this  process,  Caruso  gives  the  magnitude  of  the  intercept 

as  -Jlog  4K-.  K~  A   ,  with  K„  the  formation  constant  of  the 

96 
homocon jugate.   Utilizing  Petrella's  estimate    of  the 

limiting  conductance  for  the  ion  "pair  (IP)  H   HC1„,   of  13, 

6     -1 
K   is  calculated  to  be  4.72  x  10   mol    L ,  or  a  pK,  of  6.67. 
a  q 

The  Shedlovsky-Kay  extrapolation  here  gave  a  value  of  5.23, 
rather  different  from  that  arrived  at  from  the  Caruso  plot — 
unlike  the  situation  at  40 °C,  where  the  two  methods  gave 
fair  agreement.   This  may  be  due  to  the  difference  in 
dissociation  process. 
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For  HC1  in  pure  sulfolane,  too  many  approximations 

are  involved  to  leave  the  values  of  K  with  much  certainty 

a 

except  the  order  of  magnitude.   The  results  here  indicate 

a  pK ,  on  the  order  of  6.5.   This  is  in  very  sharp 

107 

disagreement  with  that  estimated  by  Coetzee,     pK  ,  =  14.5. 

One  must  seriously  question  several  of  the  assumptions  made 
by  Coetzee  in  arriving  at  that  conclusion,  not  the  least 
of  which  are  the  assumption  of  an  activity  for  the 
protonated  sulfolane  and  assuming  that  HC1„  has  no  effect 

whatever. 

96  97 
Petrella  and  his  coworkers   '    report  findings  for  HC1 

in  sulfolane  at  35°C  similar  to  those  found  in  this  study 

at  30 °C.   They  determined  a  pK, ,  by  the  same  method,  of  9. 

As  previously  pointed  out,  the  K   results  of  that  study  are 

consistently  higher  than  those  found  here,  and  this  is  no 

exception.   Recalculating  Petrella's  results,  however,  using 

the  values  they  specify,  gives  a  pK,  of  8.4,  more 

reasonable   in  light  of  the  value  determined  in  this  study , 

6.67,  although  the  discrepancy  is  still  perplexing. 

The  explanation  for  such  high  association  must  lie 

in  the  solvation,  or  lack  thereof,  of  the  ions. 

113  114 
Parker        has  shown  that  in  dipolar  aprotic  solvents 

such  as  sulfolane,  the  solvation  of  anions  will  be,  for  the 

most  part,  a  result  of  ion-dipole  interactions  coupled  with 

interactions  caused  by  the  mutual  polarizability  of  the 

anion  and  the  solvent  molecule.   Thus  the  observation  that 
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electrolytes  with  large  anions  are  more  solvated  makes  more 
sense.   By  the  same  token,  small  anions  such  as  chloride 
will  be  less  solvated;  indeed,  in  sulfolane,  chloride 
exists  as  a  naked,  and  hence  highly  reactive,  ion. 

The  proton  basicity  of  sulfolane  has  been  shown  to  be 

very  low,  on  the  order  of  -12.9  for  pK„H+ ,  making  it  one  of 

115 
the  most  protophobic  solvents.      Given  this  small  tendency 

to  solvate  the  proton,  coupled  with  highly  reactive  "naked" 

chloride  ions  and  sulfolane 's  ability  to  solvate  large 

anions,  the  formation  of  the  homocon jugate  complex  HCIq 

is  not  surprising.   Indeed,  it  has  been  shown  that  this 

species  forms  in  a  number  of  dipolar  aprotic  solvents,    ' 

4  117-119 
the  most  notable  being  acetonitrile. " '         The  naked 

chloride  ion  in  sulfolane  has  been  found  to  be  sufficiently 

-  120 
reactive  to  form  the  homoconjugate  AgCl„.. 

What  is,  perhaps,  most  surprising  is  that  at  40°C, 

the  conductance  data  fail  to  reveal  any  evidence  of 

homoconjugation .   The  Caruso  plot  shows  a  line  with  a  slope 

of  about  -0.5,  indicating  simple  dissociation,  and  the 

dissociation  constant  determined  via  that  method  agrees  well 

with  that  determined  by  the  Shedlovsky-Kay  extrapolation. 

With  almost  all  studies  of  sulfolane  done  at  30°C,  it  is 

most  difficult  to  interpret  this  fact. 

113 
As  Parker  has  pointed  out,     sulfolane  is  a  highly 

associated  liquid  at  30°C,  retaining  much  of  the  structure 

it  had  as  a  solid.   This  is  evidenced  by  the  high  viscosity. 

The  sharp  drop  in  the  viscosity  when  the  temperature  is 
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raised  to  40°C  is  indicative  of  a  decrease  in  the  struc- 
turing of  the  solvent.  This  probably  plays  a  role  in 
the  change  of  ionization  process  which  occurs  perhaps  by 
facilitating  the  solvation  of  the  anion.  This  being  the 
case,  it  is  conceivable  that  the  behavior  of  sulfolane  as 
a  solvent  at  30°C  is  not  typical  of  the  behavior  of  true 
liquid  sulfolane. 
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NM  Py-Water  Solvent  System 


Far  less  work  has  been  done  using  l-methyl-2- 

pyrrolidinone  (NM„Py)  as  a  solvent  than  has  been  done  with 

+  u     -,-p  i       -,  13,77,79,121-126    ,  .. 

the  sulfolane  solvent  systems,   '   '   '         and  there  are 

few  conductance  studies.    '      Since  both  sulfolane  and 

NM„Py  have  similar  physical  characteristics,  it  was  thought 

that  a  study  of  the  conductimetric  behavior  of  HC1  in 

NM„Py  and  NM„Py-water  solvents  would  be  useful  for 

comparative  purposes,  and  might  help  shed  additional  light 

on  the  properties  of  this  useful  solvent. 

The  conductance  results  for  this  portion  of  the  study 

appear  in  Table  II.   They  were  calculated  using  the  Bjerrum 

critical  distance  as  the  ion  size  parameter  in  the  full 

Pitts  equation.   The  results  for  pure  water  were  obtained 

69  70 
from  the  data  of  Stokes.        It  was  found  that,  for  pure 

NM„Py ,  the  expanded  equations  were  unable  to  provide  a 

satisfactory  fit  of  the  data,  consistent  with  the  general 

breakdown  of  these  equations  with  high  values  of  k<x;  for 

NM„Py,  at  a  concentration  of  0.01  M,  tea  =  0.45. 

For  those  solvents  rich  in  water  (25,  33,  and  50  M%) , 

it  was  noted  that,  as  the  HC1  concentration  was  increased, 

the  equivalent  conductances  also  increased  initially,  then 

leveled  out  and  decreased  in  the  normal  manner.   This 

anomalous  behavior  is  caused  by  a  reduction  in  the  ionic 

content  of  the  solutions,  and  in  the  past,  has  been 
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ascribed  to  adsorption  of  the  electrolyte  on  the  cell 

57 
walls.     This,  however,  is  a  time-dependent  process,  and 

means  that  the  measured  resistance  should  increase  with  time 

This  was  not  noted  to  be  the  case,  since  once  initial 

equilibrium  had  been  reached  after  an  addition,  no  drift  in 

measured  resistance  was  noted  over  a  span  of  72  hours. 

The  reason  for  this  strange  behavior  remains  unknown. 
The  conductance  parameters  were  generated  using  those 
experimental  points  obtained  after  the  maximum  in  conduc- 
tance was  reached.   The  initial  loss  of  small  amounts  of 
electrolyte  in  the  very  dilute  solutions  should  have  little 
impact  on  the  overall  conductance  results. 

Figure  29  is  a  plot  of  the  limiting  equivalent  conduc- 
tance, A  ,  vs.  the  mole  fraction  of  NM„Py  in  the  solvent. 
It  is  evident  that  additions  of  small  amounts  of  NM„Py  to 
water  cause  an  extremely  large  drop  in  the  limiting  conduc- 
tance of  HC1 ,  quite  similar  to  the  behavior  observed  for 
sulfolane;  the  decrease  in  conductance  due  to  NM^Py  is 
greater  still  than  that  observed  for  sulfolane. 

As  was  previously  pointed  out,  the  reason  for  such  a 
large  decrease  in  A   with  the  addition  of  the  nonaqueous 
component  lies  in  the  disruption  of  the  water  structure  and 
the  consequent  elimination  (for  all  practical  purposes)  of 

TO 

the  proton  jump  mechanism.   It  has  been  shown    that  NM„Py 
is  a  structure  breaker  with  respect  to  water,  yet  not  as 
much  of  a  structure  breaker  as  sulfolane.   Why  then  is  the 
effect  of  adding  NM„Py  to  water  greater  than  that  of  adding 
sulfolane? 
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As  shown  in  a  previous  section,  NM„Py  associates 
strongly  with  water.   Whereas  adding  sulfolane  to  water  is 
a  structure  breaking  process,  evidenced  by  the  negative 
curve  for  the  excess  viscosity,  addition  of  NM„Py  to  water 
is  a  structure  making  process.   This  can  be  seen  in  the 
large  positive  excess  viscosity.   While  the  water  structure 
is  being  disrupted,  the  over-all  structuring  in  the  solvent 
is  on  the  increase.   This  structure,  however,  evidently  does 
not  lend  itself  to  a  proton  jump  mechanism,  and  thus  the 
limiting  conductance  decreases.   In  addition,  the  large 
increase  in  viscosity,  in  and  of  itself,  helps  decrease  the 
limiting  conductance  by  decreasing  the  ionic  mobility. 
These  two  facts  together  indicate  that  the  structure  within 

the  solvent  is  probably  the  NM„Py-2H90  entity,  with  little 

13 
association  between  the  individual  entities.     Were 

aggregates  of  these  entities  to  form,  the  probability  of 

proton  jump  would  be  increased,  and  this  is  not  noted. 

Figure  30  is  a  plot  of  the  pK,  (=  log  K„ )  of  HC1  vs. 
the  mole  fraction  of  NM„Py  in  the  solvent.   The  outstanding 
feature  of  this  plot  is,  first,  its  over-all  linearity,  and 
second,  the  lack  of  any  sharp  increase  at  pure  NM„Py,  as 
was  noted  for  sulfolane.   This  is  indicative  of  the  essential 
difference  between  the  two  solvents. 

As  noted  previously,  sulfolane,  with  a  pK„„+  of  -12.9, 
is  highly  protophobic,  lending  minimal  solvation  stability 
to  the  proton.   NM„Py ,  however,  is  quite  different  in  this 
respect .   Previous  studies  have  shown  that  NM^Py  is  far  more 
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127 
basic,  with  respect  to  the  proton,  than  sulfolane,     which 

would  indicate  more  of  a  tendency  to  solvate  the  proton. 

The  measured  pKgH+  for  NM„Py  is  -0.17,  indicative  of  the 

fact  that,  unlike  the  protophobic  sulfolane,  NM„Py  can  and 

does  protonate  readily.   This  ability  of  the  solvent  to 

protonate  (and  thus  stabilize  the  proton)  is  the  principal 

reason  that,  in  NM„Py,  HC1  is  a  moderately  strong  acid, 

with  K  =  161.46  mol  '  L.   The  very  low  limiting  conductance 
a 

of  HC1  in  pure  NM„Py  is  probably  due  to  the  lack  of  a 
proton  jump  mechanism  coupled  with  the  high  viscosity,  even 
more  so  than  for  the  NM„Py-water  mixtures.   As  with 
sulfolane,  anion  solvation  is  expected  to  be  small,  although 
greater  than  in  sulfolane  due  to  a  slight  positive  charge 

+  v,     •+■         128 

on  the  nitrogen. 

The  linear  dependence  of  pK,  on  mole  fraction  NM„Py 

in  the  middle  region  of  Figure  30  can  probably  be  attributed 

to  the  greater  stabilization  (by  solvation)  of  the  proton 

with  increasing  water  content.   This  solvation  would  be 

primarily  by  the  species  NM„Py • 2H„0,  which,  apparently,  is 

well  able  to  stabilize  the  proton;  in  the  solvent  containing 

50  M%  NM„Py,  the  K   is  only  10.40  mol-1  L,  and  by  some 
^         a 

treatments,    this  is  considered  almost  completely  dis- 
sociated.  For  solvents  with  less  than  this  mole  fraction 
of  NM„Py ,  HC1  exists  primarily  as  a  strong  electrolyte. 

For  HC1  in  pure  NM„Py,  the  value  of  pK_  was  found  to 

79 
be  2.8  by  several  methods,    in  good  agreement  with  the  value 

of  2.2  obtained  in  this  study.   By  contrast,  the  results  of 
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126 
Emelin  and  his  coworkers,     yielding  a  value  for  pK   of 

a 

4.08,  are  subject  to  question.   Although  their  values  for 
the  equivalent  conductances  of  HC1  in  NM„Py  are  reasonable 
in  light  of  the  values  obtained  in  this  study,  it  is 

difficult  to  see  how  they  arrived  at  a  value  for  A   of  33 

9     _i 

S  era"  mol   .   The  value  obtained  in  this  study,  utilizing 

2     -1 

the  full  Pitts  equation,  was  8.848  S  cm  mol   ,  quite 

reasonable  in  view  of  the  constituent  equivalent  conduc- 
tances.  A  value  of  33  suggests  a  conductimetric  behavior 
of  Type  V,  and  this  is  clearly  not  the  case. 

Once  again,  the  electrostatic  theories  of  association 
show  their  weakness.   For  those  systems  where  HC1  is  a 
moderately  strong  electrolyte ,  the  Bjerrum  theory  of 
association  can  provide  a  fair  prediction  of  the  association 
constant.   Where  HC1  is  a  weak  electrolyte,  the  theory 
breaks  down.   For  0- ,  25-,  33-,  and  50-M%  NMgPy-water 
solvents,  the  Bjerrum  K  '  s  are,  respectively,  1.28,  3.31, 
4.31,  and  6.73  mol-1  L.   At  75  M%,  the  predicted  K   of 
12  mol"   L  compares  poorly  with  the  experimental  value  of 
60.63  mol    L. 

It  is  quite  interesting  to  note  that  the  method  of 

variation  of  the  ion  size  parameter,  mentioned  previously, 

was  found  to  work,  albeit  poorly,  for  those  systems  where 

the  Bjerrum  theory  provided  a  fair  approximation  of  the 

experimental  K  .   This  was  also  found  for  the  other  solvent 
^  a 

systems  studied,  as  already  mentioned. 
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There  would  thus  appear  to  be  a  positive  correlation 
between  the  ability  of  the  Bjerrum  theory  to  predict  the 
association  constant  and  a  minimum  in  a.  with  variation  of 
a.   The  Bjerrum  theory  best  predicts  the  behavior  of 

moderately  strong  electrolytes ,  and  for  these  systems 

57 
variation  of  a   does  produce  a  minimum  in  a..     This  may 

in  part  explain  some  of  the  confusion  concerning  the  value 

chosen  for  a   that  has  been  seen  in  the  past.   For  weak 

electrolytes,  where  the  Bjerrum  theory  cannot  predict  the 

experimental  K  ,  the  method  of  variation  of  a   does  not  work, 
a 

and  thus,  if  the  system  under  study  includes  a  weak 
electrolyte,  some  other  method  for  arriving  at  a  value  for 
a,  such  as  the  Bjerrum  distance  used  in  this  study,  should 
be  used. 


CHAPTER  VII 
CONCLUSIONS 

Both  sulfolane  and  NM„Py  are  dipolar  aprotic  solvents 
with  similar  physical  properties,  yet  they  differ  in  one 
major  aspect.   Sulfolane,  as  reflected  by  its  very  low  proton 
basicity,  is  protophobic,  while  NM„Py  can  hardly  be 
classified  as  such.   This  difference  is  manifested  in  the 
differing  conductimetric  behavior  of  HC1  in  the  two  solvents. 

The  limiting  conductance  of  HC1  was  found  to  drop 
drastically  with  increasing  sulfolane  content  in  the 
solvent.   This  is  attributed  to  the  disruption,  by  sulfolane, 
of  the  water  structure  that  supports  the  proton  jump 
mechanism  responsible  for  the  high  conductance  of  strong 
acids  in  water.   The  strength  of  HC1 ,  as  measured  by  the 
dissociation  constant,  was  also  found  to  decrease  with 
increasing  sulfolane  content,  but  in  a  more  moderate  manner. 
The  exception  to  this  was  in  pure  sulfolane,  where  the  HC1 
was  found  to  exist  primarily  in  an  associated  form. 

The  full  Pitts  equation,  using  the  Bjerrum  critical 
distance  as  the  ion  size  parameter,  was  able  to  provide  a 
satisfactory  fit  of  the  conductance  data  for  HC1  in  the 
sulf olane-water  solvents  at  all  temperatures  studied.   The 
exception  to  this  was  HC1  in  pure  sulfolane,  at  30  and  40°C. 
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These  data  were  analyzed  using  the  method  of  Caruso.   At 
30°C,  homoconjugation  was  found,  and  prior  work  by  Benoit 
has  shown  that  the  species  formed  is  HClI.   At  40°C,  no 
evidence  of  homoconjugation  was  observed,  and  the  conduc- 
tance data  were  analyzed  by  using  the  extrapolation  method 
of  Shedlovsky  and  Kay.   It  is  possible  that  this  change  in 
the  dissociation  process  with  temperature  is  caused  by  a 
decrease  in  the  structuring  of  the  sulfolane,  which  at  30°C 
retains  much  of  the  structure  it  had  as  a  solid.   The 
homoconjugation  at  30°C  is  probably  caused  by  the  ability 
of  sulfolane  to  solvate  large  ions  coupled  with  very  poor 
solvation  of  the  anion,  CI  .   It  is  known  that  a  number  of 
homoconjugates  such  as  AgCl„  form  in  sulfolane. 

Since  almost  all  previous  work  with  sulfolane  has  been 
done  at  30°C,  it  is  difficult  to  know  what  other  effects 
(if  any)  are  observed  with  rising  temperature.   The  results 
of  this  study  suggest  that  at  30°C  sulfolane  may  not  be  a 
typical  liquid  dipolar  aprotic  solvent,  and  further  work 
at  elevated  temperatures  is  needed. 

The  conductimetric  behavior  of  HC1  in  NM9Py-water 
solvents  rich  in  water  is  very  similar  to  the  behavior  noted 
in  sulfolane-water  mixtures.   The  limiting  conductance  of 
HC1  fell  rapidly  with  increasing  NM„Py  content,  attributable 
to  the  disruption  of  the  water  structure  by  NM~Py .   Both 
NMpPy  and  sulfolane  have  been  found  to  be  water  structure 
breakers,  but  it  was  noted  that,  whereas  sulfolane  caused  a 
general  decrease  in  solvent  structuring,  NM^Py  caused 
a  marked  increase  in  solvent  structuring,  attributable 
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to  the  interaction  between  NM9Py  and  water  to  form  the 
associated  species  NM„Py • 2H„0. 

The  ion  association  present  in  HC1  solutions  rose 
gradually  with  increasing  NM„Py  content  in  the  solvent,  but 
even  at  50  M%  NM„Py  could  be  considered  only  slightly 
associated.   HC1  was  found  to  exist  as  a  moderately  strong 
acid  in  pure  NM„Py .   This  may  be  attributed  to  the  ability  of 
NM„Py  to  enter  into  hydrogen  bonding  at  the  carbonyl  oxygen , 
and  thus  stabilize  the  proton  by  solvation.   The  low 
limiting  conductance  found  for  HC1  in  NM„Py  and  NM„Py-water 
solvents  is  believed  to  be  caused  by  the  lack  of  a  proton 
jump  mechanism,  coupled  with  a  high  solution  viscosity, 
and  the  possible  formation  of  discrete  hydrodynamic  entities. 
As  for  sulfolane,  and  most  other  dipolar  aprotic  solvents, 
solvation  of  the  anion  is  expected  to  be  poor. 

The  80  W%  MetOH-H20  solvent  was  found  to  foster 
slight  association  in  the  electrolytes  studied,  which  are 
usually  considered  to  be  strong  electrolytes.   NaCl  and  HC1 
were  found  to  be  slightly  associated,  while  NaAc  and 
Tris-HCl  were  found  to  be  associated  to  a  slightly  greater 
degree.   HAc  was  found  to  be  strongly  associated,  with  a 
dissociation  constant  reasonable  in  view  of  that  found  in 
similar  solvent  mixtures  using  other  alcohols.   The  larger 
pK,  found  for  the  80  W%  MetOH-H„0  solvent  may  be  caused  by 
steric  hindrance  in  the  solvation  of  the  proton,  resulting 
from  the  larger  size  of  the  methoxyethanol  molecule.   The 
proton  jump  mechanism  is  noted,  but  in  a  very  reduced  form, 
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and  this  has  been  explored  by  Thun . 

Tris-HCl  was  also  studied  in  water.   It  was  found  to 
be  a  strong  electrolyte,  essentially  dissociated.   The 
limiting  conductance  for  the  TrisH   ion  was  found  to  be 
smaller  than  expected.   This  is  consistent  with  the  ability 
of  the  TrisH   ion,  with  three  OH  groups,  to  hydrogen  bond 
with  the  water  structure.   Transfer  to  the  80  W%  MetOH-H„0 
solvent  caused  an  increase  in  ion  association  and  a  decrease 
in  the  limiting  conductance.   This  is  reasonable  in  view  of 
the  lesser  ability  of  the  mixed  solvent  to  solvate  the 
TrisH   ion. 

It  was  noted  that,  for  those  systems  where  the  electro- 
static theory  of  ion  association  of  Bjerrum  was  able  to 
provide  a  fair  approximation  of  the  association  constant, 
a  minimum  in  the  standard  deviation  of  the  conductance  could 
be  obtained  by  varying  the  ion  size  parameter.   For  those 
systems  where  the  Bjerrum  theory  failed,  so  did  the  variation 
method.   As  a  result,  to  be  consistent,  the  Bjerrum  critical 
distance  was  used  throughout  this  study.   Where  the  systems 
studied  have  moderate  association  approximated  by  the  Bjerrum 
theory,  the  variation  method  of  treatment  appears  to  have 
merit . 


APPENDICES 

APPENDIX  A CONDUCTANCE  DATA 

APPENDIX  B COMPUTER  PROGRAMS 


APPENDIX  A 
CONDUCTANCE  DATA 


SOLVENT 25  tt%    Sulfolane-Water 

TEMPERATURE 25   C 

SOLUTE HC1 


c(mol  L  ) 
3.4981E-04 
6.2844E-04 
9.9713E-04 
1.4132E-03 
1.8148E-03 
2.2121E-03 
2.6116E-03 
3-3772E-03 
3.7757E-03 
4.1563E-03 
4.5339E-03 


A(S  cm2 

mol   ) 

103. 

■  099 

102. 

,649 

102. 

.179 

101, 

•  979 

101. 

,470 

101. 

,180 

100. 

,900 

100. 

,420 

100. 

,210 

99. 

■  921 

99. 

771 

SOLVENT 25  M%    Sulfolane-Water 

TEMPERATURE 30   C 

SOLUTE HC1 


c(mol  L  ) 
1.38H6E-03 
2.4728E-03 
3-7505E-03 
4.8106E-03 
6.2202E-03 
7.7413E-03 
8.8779E-03 
1.0112E-02 
1.1308E-02 
1.2536E-02 
1.4135E-02 
1.5442E-02 


A(S  cm" 

mol   ) 

112, 

.469 

111, 

•  707 

110, 

•  919 

110, 

.347 

109. 

,577 

108. 

,859 

108. 

,382 

107. 

,881 

107- 

,454 

106. 

,980 

106. 

447 

106. 

092 
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SOLVENT 25    MJI    Si 

RATURE 40 

SOLUTE HCl 


TEMPERATURE 40  °C 


c(mol  L   ) 

2.3435E-03 
3.8634E-03 
5-9749E-03 
7.3956E-03 
9.4757E-03 
1.1179E-02 
1.2948E-02 
1.4602E-02 
1.6615E-02 
1.9484E-02 
2.2091E-02 
2.5953E-02 


A(S  cm   mol   ) 

131. 

■  493 

130. 

,056 

128. 

.391 

127. 

,477 

126. 

,132 

125. 

,114 

123- 

■  990 

123- 

,130 

122. 

,140 

121. 

,128 

120, 

,000 

118. 

,656 

SOLVENT 50  l\%    Sulfolane-Water 

TEMPERATURE 25  °C 


SOLU1 

'E HCl 

c(mol  L   ) 

2     -1 
A(S  cm   mol   ) 

2.4543E-03 

33-647 

2.7605E-O3 

33.489 

3.1395E-03 

33.234 

3.5024E-03 

33-037 

4.2715E-03 

32.608 

5.1971E-03 

32.153 

5.7115E-03 

31.899 

6.6713E-03 

31.446 

7.1798E-03 

31.255 

7.7786E-03 

30.969 

-152- 
VENT— 

TURK- 
SOLUTE HC1 


SOLVENT 50  112  Sulfolane-Water 

TEMPERATURE 30   C 


c(mol  L  ) 
1.7056E-03 
3.1260E-03 
4.6395E-03 
6.3705E-03 
7.8767E-03 
9.33721E-03 
1.0844E-02 
1.2363E-02 
1.4034E-02 
1.5823E-02 
1.7826E-02 
2.0621E-Q2 


A(S  cm 

2   i-^ 
mol   ) 

39. 

,130 

37- 

,994 

36. 

■  750 

35. 

,624 

34. 

,770 

33. 

•  994 

33. 

•  385 

32, 

•  747 

32. 

,136 

31. 

■  507 

30. 

,889 

30. 

,130 

SOLVENT 50  M   Sulfolane-Water 

TEMPERATURE 40   C 

SOLUTE HC1 


c 

(mol  L   ) 

1. 

■8125E-03 

2. 

.0720E-03 

2. 

•5791E-03 

3. 

,  1078E-03 

3. 

,  5393E-03 

4. 

.0950E-03 

4. 

■5835E-03 

5. 

•0959E-03 

5. 

,  5862E-03 

1 

A(S  cm' 

mol   ) 

46, 

,188 

45. 

,894 

45- 

,296 

44. 

,809 

44. 

•  376 

43. 

,849 

43- 

,423 

43. 

,022 

42. 

,772 
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SOLVENT— 75  n%   Sulfolane-Water 

RATURS 25 

SOLUTE HC1 


TEMPERATURE 25  °C 


c(mol  L   ) 

1.3943E-03 
2.6590E-03 
4.1185E-03 
5.8971E-03 
7-8847E-03 
1.0373E-02 
1.2904E-02 
1.5524E-02 
1.8919E-02 
2.1870E-02 
2.4834E-02 


A(S 

2            -1 
cm     mol      ) 

14. 

.060 

12, 

,644 

11, 

.314 

10. 

.172 

9. 

.284 

8. 

,461 

7- 

,810 

7. 

305 

6. 

,768 

6. 

436 

6. 

361 

SOLVENT 75  Wk    Sulfolane-Water 

TEMPERATURE 30   C 

SOLUTE HC1 


c(mol  L   ) 

1.0817E-03 
1.8025E-03 
2.5500E-03 
3.4895E-03 
4.6936E-03 
5.9074E-03 
7.2180E-03 
8.8548E-03 
1.0214E-02 
1.1791E-02 
1.3485E-02 
1.5766E-02 


c 

A(S    cm' 

mol      ) 

15. 

,607 

14. 

,830 

13. 

,694 

12. 

,676 

11. 

,747 

10. 

,980 

10. 

,301 

9- 

,669 

9. 

,194 

8, 

.758 

8. 

.343 

7. 

.872 
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SOLVENT 75  Mi  Sulfolane-Water 

TEMPERATURE 40   C 

SOLUTE HC1 


c(mol  L   ) 

1.3318E-04 

4.0819E-04 

7.6171E-04 

1.1057E-03 

1.4410E-03 

1.8107E-03 

2.2073E-03 

2.5469E-03 

2.8866E-03 

3.23HE-03 

3.5768E-03 

3.9077E-03 

4.2309E-03 

4.5565E-03 

4.9024E-03 


2 

:   ,-1, 

A(S  cm 

mol   ) 

25- 

324 

23- 

608 

22. 

133 

21. 

039 

20. 

087 

19. 

069 

18. 

,214 

17. 

,627 

17. 

,123 

16, 

.612 

16, 

.166 

15 

.788 

15 

.430 

15 

.108 

14 

.766 

SOLVENT 85  M*  Sulfolane-Water 

TEMPERATURE 25   C 

SOLUTE KC1 


c(mol  L   ) 

2.  1200E-03 
4.6748E-03 
9.7251E-03 
1. 1215E-02 
1.2971E-02 
1.4923E-02 
1.7203E-02 
1.8653E-02 
2.0380E-02 

2.2615E-02 
2.  7H3E-02 
3-3252E-02 


A(S  cm 

2     -1 

mol   ) 

6, 

.822 

5, 

.173 

3, 

.981 

3, 

•  723 

3, 

,506 

3, 

.318 

3. 

,165 

3. 

,059 

2. 

,955 

2. 

,838 

2. 

,655 

2. 

,461 
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S0LVEI1T 85    n%   Sulfolane-Water 

TEMPERATURE 30   C 


SOLUTE HG1 


c(mol  L  X)  A(5  cm2  mol  -1 ) 

2.3228E-03 
3.9677E-03 
6.6773E-03 
8.9155E-03 
1.1130E-02 
1.5784E-02 
1.8941E-02 
2.3524E-02 
2.4246E-02 
2.8467E-02 
3.4505E-02 
5.2333E-02 


7- 

■  633 

6. 

■  503 

5. 

•  133 

4, 

■  575 

4. 

,218 

3. 

,663 

3. 

.541 

3. 

,176 

3. 

,111 

3. 

,096 

2. 

■  790 

2. 

,448 

SOLVENT 85  M%   Sulf olane-Water 

TEMPERATURE >i0      C 

SOLUTE HC1 


c(mol  L  ) 
2.3419E-03 
3.3141E-03 
4.5002E-03 
5.4610E-03 
6.5H2E-03 
7.6747E-03 
8.5836E-03 
9.7424E-03 
1.1434E-02 
1.2571E-02 
1.4119E-02 


A(S 

cm   mol   ) 

8. 

,741 

7. 

■  572 

6, 

,761 

6. 

,269 

5. 

,872 

5. 

.439 

5. 

,197 

H. 

■  939 

4. 

,676 

h. 

•  515 

'4. 

,340 
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SOLVENT 100  M%    Sulfolane 

TEMPERATURE 30   C 

SOLUTE HC1 


c(mol  L  ) 
1.7803E-03 
4.5526E-03 
5.4460E-03 
6.8664E-03 
8.4110E-03 
1.1232E-02 
1.4040E-02 
1.7210E-02 
1.9644E-02 
2.2716E-02 


A(S 

2           n"^ 
cm      mol       ) 

0. 

,267 

0. 

,238 

0. 

,234 

0, 

,225 

0. 

,222 

0. 

,223 

0, 

,224 

0. 

,226 

0. 

,227 

0. 

,228 

SOLVENT 100  K%    Sulfolane 

TEMPERATURE 40  °C 

SOLUTE HC1 


c(mol  L  ) 
5.4732E-04 
1.0673E-03 
1.6745E-03 
2.2722E-03 
3.2895E-03 
4.5549E-03 
5.6198E-03 
6.6298E-03 
7.7673E-03 
9.6316E-03 


A(S 

2         .-1. 

cm      mol       ) 

0. 

.812 

0. 

■  511 

0. 

,384 

0. 

,308 

0. 

,239 

0. 

,196 

0. 

,182 

0. 

,163 

0. 

,153 

0. 

,136 
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SOLVENT 25  m    NM2Py-Water 

TEMPERATURE 25   C 

SOLUTE HC1 


c(mol  L  ) 
2.042HE-03 
4.2705E-03 
5.4351E-03 
7.1128E-03 
8.2325E-03 
9-3255E-03 
1.0622E-02 
1.1926E-02 
1.3467E-02 
1.5137E-02 
1.7001E-02 
1.8770E-02 


A(S    cm' 

mol       ) 

49. 

.751 

50. 

,230 

50. 

.275 

50. 

,047 

49. 

,886 

49. 

,403 

49. 

,185 

48, 

•  990 

48, 

,757 

48, 

,512 

48, 

.243 

48, 

,040 

SOLVENT 33  M$  NM2Py-Water 

TEMPERATURE 25   C 

SOLUTE HC1 


c(mol  L  ) 
1.H722E-03 
3.3547E-03 
4.8889E-03 
6.4344E-03 
8.0910E-03 
9-9425E-03 
1.1353E-02 
1.2801E-02 
1.4247E-02 
1.5845E-02 
1.7423E-02 


A(S 

t 

cm' 

mol       ) 

32, 

,784 

34. 

,744 

35. 

,259 

35. 

•  320 

35. 

,312 

35- 

,106 

35. 

,004 

34. 

,884 

34. 

,773 

34. 

,647 

34. 

■  513 
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SOLVENT 50    m    NM2Py-Watf 

TEMPERATURE 25   C 

SOLUTE HC1 


c(mol  L  ) 
6.2246E-04 
1.0941E-03 
1.9876E-03 
3.0661E-03 
4.0665E-03 
5.3296E-03 
6.6524E-03 
7-9587E-03 
9.6281E-03 
1.113SE-02 
1.2467E-02 
1.3784E-02 


A(S    cm' 

mol       ) 

20. 

.428 

22. 

,514 

23. 

•  507 

23. 

,628 

2°.. 

.492 

23. 

.242 

22, 

.945 

22, 

.654 

22, 

.218 

21, 

.920 

21, 

.675 

21, 

.441 

SOLVENT 75  M#  NM2Py-Water 

TEMPERATURE 25   C 

SOLUTE HC1 


c(mol  L   ) 

1.0018E-03 
1.8847E-03 
2.9497E-03 
3.9830E-03 
4.8933E-03 
5.8259E-03 
6.6806E-03 
7.6254E-03 
8.5401E-03 
9.5531E-03 


A(S 

1 

cm' 

2           -1 
mol       ) 

14, 

.397 

13, 

.346 

12, 

.762 

12, 

.342 

11, 

.934 

11, 

.627 

11, 

.287 

10, 

■  998 

10, 

,744 

10, 

.448 

-159- 


SOLVENT 100  m    NM2Py 

TEMPERATURE 25   C 

SOLUTE HC1 


c(mol  L  ) 
1.9017E-03 
3.0819E-03 
4.0778E-03 
5.1046E-03 
6.1461E-03 
7.0249S-03 
7.9772E-03 
9.0077E-03 
1.0449E-02 
1.1856E-02 
1.3139E-02 
1.4758E-02 


A(S 

cm 

2     -1 
mol   ) 

7. 

003 

5. 

558 

4. 

■  922 

4. 

■  493 

4. 

,190 

3. 

■  995 

3. 

,827 

3. 

,681 

3. 

■  521 

3. 

.399 

3. 

,301 

3. 

,211 

SOLVENT 80  Vl%   MetOH-Water 

TEMPERATURE 25   C 


SOLUTE HC1 


c(mol  L   ) 

1.03O8E-O3 

2.3159E-03 
3.6672E-03 
4.9493E-03 
6.8956E-03 
9-5516E-03 
1.1168E-02 
1.2813E-02 
1.4667E-02 
1.6899E-02 
1.8978E-02 
2.2214E-02 


A(S  cm' 

mol   ) 

50. 

.235 

48. 

,726 

47- 

,695 

46, 

,880 

45. 

,898 

44. 

,851 

44. 

,328 

43. 

,833 

43. 

,337 

42, 

,881 

42. 

,455 

41, 

,892 
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SOLVENT 80  W2  MetOH-Water 

te: 

V1PERATURE 25   C 

SOLUTE NaCl 

c(mol  L   ) 

A(S  cm2  mol-1) 

8.8831E-04 

29.752 

1.6551E-03 

29.035 

2.3491E-03 

28.5^2 

3.1126E-03 

28.098 

3.8850E-03 

27.727 

4.7360E-03 

27-357 

5.7255E-03 

26.999 

6.6050E-03 

26.707 

7.5226E-03 

26.433 

8.6949E-03 

26.117 

9.9230E-03 

25.823 

1.1851E-02 

25.418 

SOLVENT 80  MetOH-Water 

TEMPERATURE 25   C 

SOLUTE NaAc 


c(mol  L  ) 
9.3058E-04 
1.2630E-03 
1.6773E-03 

2.0625E-03 
2.6011IE-03 
3.7400E-03 
4.6878E-03 


5.6674E-03 
6.6449E-03 
7.3989E-03 
8.4268E-03 

9.8518E-03 


A(S  cm 

2     -1 
mol   ) 

23. 

■  958 

23. 

,676 

23. 

.362 

23. 

,101 

22, 

,773 

22, 

,183 

21, 

,768 

21, 

•  397 

21, 

,069 

20. 

,839 

20, 

,554 

20, 

,194 

-161- 

SOLVEMT 80    \n   MetOH-Water 

TEMPERATURE 25   C 

SOLUTE TrisHCl 


-1  2     -1 
c(mol  L   )                  A(S  cm   mol   ) 

1.0193E-03  24.868 

2.4322E-03  23-686 

3.7122E-03  22.927 

4.5503E-03  22.523 

5.4897E-03  22.126 

6.2258E-03  21.845 

7.0050E-03  21.576 

8.1082E-03  21.235 

9.3267E-03  20.892 

1.0713E-02  20.541 

1.1796E-02  20.293 

1.3009E-02  20.038 


SOLVENT 80  W2  MetOH-Water 

TEMPERATURE 25   C 

SOLUTE HAc 


-1  2     -1 
c(mol  L   )                  A(S  cm   mol   ) 

1.4286E-03  0.033 

2.5720E-03  0.033 

3.9471E-03  0.035 

5-3867E-03  0.032 

6.6175E-03  0.032 

8.1364E-03  0.031 

9.2652E-03  0.030 

1.1008E-02  0.029 

1.2341E-02  0.029 

1.3441E-02  0.029 
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SOLVENT- 

— Water 

TEMPERATURE  - 

—  25  °C 

SOLUTE- 

— Tris-HCl 

c(mol  L   ) 

A(S  cm2  mol  -1) 

3.8086E-04 

104.878 

4.3331E-04 

104.808 

6.4097E-04 

103.779 

9.8939E-04 

103.469 

1.2467E-03 

102.929 

1.4537E-03 

102.809 

1.5429E-03 

102.809 

1.8155E-03 

102.869 

1.9323E-03 

102.239 

2.1055E-03 

101.810 

2.2491E-03 

101.760 

2.2803E-03 

101.940 

2.3112E-03 

102.279 

2.9637E-03 

101.590 

2.9674E-03 

101.310 

3.4152E-03 

100.700 

3.6500E-03 

100.620 

3.7873E-03 

100.890 

4.8032E-03 

99.789 

4.8078E-03 

99.883 

6.3979E-03 

99. 032 

6.6171E-03 

98.838 

7.4038E-03 

98.502 

7.7206E-03 

98.449 

8.3647E-03 

98.095 

8.7773E-03 

97.997 

9.3236E-03 

97.721 

9.5604E-03 

97.712 

9.9768E-03 

97.467 

1.0202E-02 

97.449 

1.1058E-02 

97.066 

1.1061E-02 

97.139 

1.1907E-02 

96.798 

1.2002E-02 

96.831 

1.2711E-02 

96.902 

1.2879E-02 

96.560 

1.3707E-02 

96.237 

1.3735E-02 

96.292 

1.4884E-02 

95.993 

1.4908E-02 

95.909 

1.5952E-02 

95.615 

1.6109E-02 

95.633 
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SOLVENT- 

— Water 

TEMPERATURE- 

—  37  °C 

SOLUTE- 

— Tris-HCl 

c(mol  L   ) 

A(S  cm2  mol  X) 

7.1244E-04 

131.035 

1.00S1E-03 

130.355 

2.1840E-03 

128.506 

2.6237E-03 

127.947 

3.7574E-03 

126.977 

4.6154E-03 

126.257 

5.0611E-03 

125.968 

6.1003E-03 

125.288 

7.2061E-03 

124.628 

7.8365E-03 

124.308 

8.8483E-03 

123.689 

9.4913E-03 

123.399 

1.0604E-02 

122.859 

1.0870E-02 

122.739 

1.1411E-02 

122.509 

1.1810E-02 

122.319 

1.2318E-02 

122.090 

1.2846E-02 

121.890 

1.3355E-02 

121.740 

1.3761E-02 

121.590 

1.4954E-02 

121.090 

1.4990E-02 

121.110 

1.5620E-02 

120.850 

1.5660E-02 

120.890 

APPENDIX  B 
COMPUTER  PROGRAMS 


PITTSEQM/BAS 


10  'PROGRAM  3Y  S.  DAVID 

FLA.,  GAINESVILLE,  FL 
20  CLS 
30  PRINT"  *** 

40  PRINT :PRINT"THIS  IS 

INTER  ONLY!" 

50  PRINT: PRINT"BE  SURE 

RECISION!" 

60  PRINT: PRINT" COUNTER 

MENT  1090;  COUNTER  'M' 

N  STATEMENT  10 80 

ONS  IN  STATEMENT 


KLEIN 
32611 


DEPT 
USA 


OF  CHEM. 


UNIV.  OF 


THIS  IS  THE  PITTS  EQUATION  PROGR 

A  REVISION.  RESULTS  OUTPUT  TO  PR 

TO  ENTER  ALL  NUMBERS  IN  DOUBLE  P 

!J'  SETS  #  OF  DATA  SETS  IN  STATE 
SETS  #    OF  ION  SIZE  ITERATIONS  I 
COUNTER  'K1    SETS#  OF  LAM-0  ITERATI 
1060.": PRINT 


70  DEFDBL  A-H,0,Q-Z 

80  DIM  EC(25),EL(25),XL(25),QL(25),AL(25)>CI(25),RC(25) 

,YP(25),YQ(25),YR(25)5EP(25),EQ(25),ER(25),TP(25),SP(25 

),RL(25),D(25),Y(25),G(25),FA(25)JQK(25)JX(25),BL(25)5C 

L(25),DF(25),DS(25),TA(25),TB(25),TC(25),AK(25) 

90  J=0 

100  J=J+1 

110  INPUT"INPUT  DC, TEMP(K), VISC. , DENS. ,MW  SOLV. ,MW. SOLU 

.";DC,T,V,DE,MS#,MV# 

120  PRINT" NOW  INPUT  THE  #  OF  POINTS, THE  APPROX.  LAMBDA- 

0" 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 


VIS; 


";V 


MWSOLU=";MV# 


EL=####.# 


INPUT  N,OL 

PRINT" INPUT  THE  CONCENTRATION, CONDUCTIVITY  PAIRS 

FOR  1=1  TO  N 

INPUT  EC(I),EL(I) 

EL(I)=0.999505DO»SL(I) 

NEXT  I 

LPRINT"N=";N;";  DC=";DC;";  T(K)=": 

LPRINT"DEN3=" ;DE; " ;  MWSOLV=" ;MS# ; ' 

LPRINT 

FOR  1=1  TO  N:LPRINTUSING"EC=##. ####[[[[ 
###" ;EC(I) ,EL(I) :NEXT:LPRINT 
230  DB=DC*T 
240  Z=DB:GOSUB  1600 
250  A1=3.643252D6/(DB^Z3) 
260  PRINT  "ENTER  AO  FOR  THIS  RUN.'" 
270  INPUT  AO 

280  AP=1.158266D6/(DB*Z3) 
290  BP=2.345787D13/(DB*DB*DB) 

300  GP=82.34l64DO/(V*(Z3)) 

310  HP=2.796299D6/(DB*Z3) 

320  PRINT"INPUT  THE  AO  INCREMENT."  :  INPUT  AI 

330  M=0 

340  Z9=Z3 

350  M=M+1 

360  A0=A0+AI 

370  LPRINT" A0=" ;A0 
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380  EA=5-O2O330D9*AO/(Z9) 

390  Z=2.0D0:GOSUBl600:RT=Z3 

400  PRINT" AP=" ; AP , "3P= " ; BP , " GP=" ; GP , "HP=" ; HP , " RT= " ; RT , " 

BA=";BA 

410  PRINT: PRINT 

420  K=0 

430  K=K+1 

440  FOR  1=1  TO  N 

450  XL(I)=OL 

460  QL(I)=XL(I) 

470  AL(I)=EL(I)/QL(I)  :  CI( I) =AL(I)*EC(I)  :  Z=CI(I):GOS 

UB  1600:RC(I)=Z3 

480  YP(I)=BA*RC(I)  :  YQ( I ) =1 . 707107D0*YP( I)  :  YR(I)=2.7 

07107D0*YP(I) 

490  Q=YP(I) :GOSUB  1100 

500  EP(I)=RE 

510  Q=YQ(I)  :  GOSUB  1100 

520  EQ(I)=RE 

530  Q=YR(I)  :  GOSUB  1100 

540  ER(I)=RE 

550  Z=YP(I) :GOSUB  1500 :TA ( I) =Z2 

560  Z=YQ(I) : GOSUB  1500 :TB( I) =12 

570  Z=YR(I) :GOSUB  1500 :TC (I) =Z2 

580  TP(I)=0.75DO/(RT+YP(I))+(3.0DO/8.0DO)*(3.0DO*RT*TB( 

I)*EQ(I)/(RT+YP(I))-2.0D0*TA(I)*EP(I)) 

590  ZK=((1.0D0+YP(I))*(RT+YP(I))) 

600  SP(I)=(9D0*RT-10D0+YP(I)*(3D05RT+1D0)+2D0*(YP(I)*YP 

(I)))/(8D0*(ZK*ZK))-TA(I)*EP(I)/(4D0*(1D0+YP(I)))+(7D0* 

RT*TC(I)*ER(I))/(16D0*((1D0+YP(I))*(1D0+YP(I)))*(RT+YP( 

I)))+(RT*TB(I)*EQ(I))/(16D0*(1D0+YP(I))*(RT+YP(I))) 

610  RL(I)=-OL*(APsRC(I)/((1.0D0+YP(I) ) * (RT+YP( I) ) ) +BP*C 

I ( I ) *  SP ( I ) ) -GP*RC ( I ) * ( 1 . 0D0/ ( 1 . 0 DO+YP ( I ) ) -AP*  RC ( I ) / ( ( ( 1 

.0D0+YP(I))*(1.0D0+YP(I)))*(RT+YP(I)))-HP*RC(I)*TP(I)/( 

1.0D0+YP(I))) 

620  XL(I)=OL+RL(I) 

630  D(I)=ABS(XL(I)-QL(I)) 

640  IP  D(I)  >  0.001  GOTO  460 

650  Y(I)=EL(I)-RL(I) 

660  ZT=((-A1*RC(I)/(1.0D0+BA*RC(I)))) 

670  Z=1.0D1:GOSUB  1390 : ZM=ZT*Z2' 

680  Z=ZM:GOSUB  1500:G(I)=Z2 

690  FA(I)=(DE/(DE+1.0D-3*EC(I)*((2.0D0*MS#)-MV#)))*(DE/ 

(DE+1.0D-3*EC(I)»((2.0DO*MS#)-MV#))) 

700  G(I)=G(I)*FA(I) 

710  QK(I)  =  (1.0DO-AL(I))./(((AL(I))*(AL(I)))*EC(I)*G(I)) 

720  X(I)=((EL(I)*EL(I))*EC(I)*G(I))/XL(I) 
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730  PRINT"I=" ;I:PRINT:PRINT"QL=";QL(I) ,"AL=";AL(I) ,"YP  = 
";YP(I),"YQ=";YQ(I),"YR=";YR(I),"EP=";EP(I),»EQ=";EQ(I) 

,"ER=";ER(I),"TP=";TP(I),"SP=";SP(I),"RL=";RL(I),"XL="; 

XL(I),"Y=";Y(I)/'G=,';G(I),"FA=";FA(I),"QK  =  ";QK(I),"X  =  "; 

X(I) :PRINT 

740  NEXT  I 

750  SX=0  :  SY=0  :  ZX=0  :  ZY=0  :  ZZ  =  0 

760  FOR  L=l  TO  N 

770  SX=SX+X(L)  :  SY=SY+Y(L)  :  ZZ=ZZ+X(L) *Y(L)  :  ZX=ZX+( 

X(L)*X(L))  :  ZY=ZY+(Y(L)*Y(L) ) 

780  NEXT  L 

79  0  PRINT"SX=";SX,"SY=";SY,"ZZ=";ZZJ"ZX=";ZX3"ZY=";ZY 

800  AB=(N*ZX)-(SX*SX) 

810  OL=((ZX*SY)-(SX*ZZ))/AB 

820  SL=((N*ZZ)-(SX*SY))/AB 

830  EZ=-SL 

840  F=((ZZ-(SX*SY)/N)*(ZZ-(SX*SY)/N))/(ZX-(.(SX)*(SX))/N 

) 

850  Z=((ZY-(( (SY)*(SY) )/N)-F)/(N-2.0D0)) :GOSUB  l600:SO= 

Z3 

860  Z=ZX/AB:GOSUB  1600 : SB=SO*Z3 

870  Z=N/AB:GOSUB  1600 : SR=SO*Z3 

880  SS=SR 

890  FOR  1=1  TO  N 

900  BL(I)=OL-OL*(AP*RC(I)/((1.0D0+YP(I) ) * ( RT+YP( I ) ) ) +BP 

*CI(I)*SP(I))-GP*RC(I)*(l.0DO/(l.OD0+YP(I) ) -AP*RC(I)/( ( 

( 1 . ODO+YP ( I ) ) * ( 1 . ODO+YP ( I ) ) ) * ( RT+ YP ( I ) ) ) -HP*RC ( I ) *TP ( I ) 

/(1.0D0+YP(I))) 

910  Z=1.0D0+4.0D0*EC(I)*G(I)*EZ:GOSUB  1600 : AK( I) =Z3 

920  AL(I)=ABS((AK(I)-1.0D0)/(2.0D0*EC(I)*G(I)*EZ)) 

930  CL(I)=AL(I)*BL(I) 

940  DF(I)=EL(I)-CL(I) 

950  DS(I)=DF(I)*DF(I) 

960  NEXT  I 

970  SD=0 

980  FOR  L=l  TO  M 

990  SD=SD+DS(L) 

1000  NEXT  L 

1010  Z=(SD/(N-1.0D0)) :GOSUB  l600:SE=Z3 

1020  LPRINT  USIHG"LAM-0=###.  ###'■*  ;  OL  ;  :  LPRINTUSING"      3 

D  LAM-0=     ##.####[[[[";SB; : LPRINTUSING"      KA=##.### 

#[[[[";EZ 

1030  LPRINTUSING"SD  KA=##. ####[[[[" ;SS ; : LPRINTUSING" 

SD  LST  SQ  LINE=##. ####[[[["; SO 

1040  LPRINT 

1050  IF  K  <  5  GOTO  430 

1060  LPRINT  "SD=" ;SE: LPRINT 

1070  IF  U   <20  GOTO350 

1080  IF  J  <  1  GOTO  100 

1090  END 
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1100  REM  *  THIS  IS  FORT  EXPI  SUBROUTINE  * 
1110  IF  (Q-i)<0  THEN  GOTO  1200 

1120  IFQ=0  THEN  GOTO  1180 

1130  R=l/Q 

1140  AU=1.0D0-R*(((R+3.377358D0)*R+2.0  252l6D0)*R+2.7094 

79D-l)/((((R*1.072553D0+5.7l6943D0)*R+6.945239D0)*R+2.5 

93888D0)*R+2.709^96D-1) 

1150  Z=-Q:G0SUE  1500 

1160  RE=AU*R*Z2 

1170  RETURN 

1180  PRINT  "SUBROUTINE  BOMB...X=0" 

1190  STOP 

1200  IF  (Q+3)<=0  THEN  GOTO  1250 

1210  AU=(( (((((7.122452D-7*Q-1.7663l5D-6)*Q+2.928l33D-5 

)*Q-2.335379D-4)*Q+1.664l56D-3)*Q-1.04l576D-2)*Q+5.5556 

82D-2)*Q-2.500000D-l)*Q+9-999999E-l 

1220  Z=ABS(Q) :GOSUB  1390 

1230  RE=Q*AU-Z2-5-772157D-l 

12  40  RETURN 

1250  IF  (Q+9)<=0  THEN  GOTO  1280 

1260  AU=l-((((5.176245D-2*Q+3.0ol037D0)*Q+3.243665Dl)*Q 

+2.244234D2)*Q+2.486697D2)/((((Q+3-995l6lD0)*Q+3.893944 

Dl)*Q+2.2638l8Dl)*Q+1.807837D2) 

1270    GOTO    1310 

1280   R=9/Q 

12  90  AU=1-R* (( (R+7. 659 824D-l)*R-7. 27 10 15 D-1)*R- 1.080693 

D0)/( (((R*2.5l8750D0+1.22927Dl)*R+5.92l405D0)*R-8. 66670 

2D0)*R-9.724216D0) 

1300  Z=-Q:GOSUB  1500 

1310  RE=AU*Z2/Q 

1320  RETURN 

1330  END 

1340  -DOUBLE  PRECISION  ROUTINES 

1350  'COPYRIGHT  TANDY  CORP.  1978 

1360  ' 

1370  * ***NATURAL  LOG-** 

1380  * 

1390  Z2=LOG(Z) :I0=0:I2=SGN(Z2) : IF  I2<0  Z=l/Z 

1400  IF  Z<1.065  THEN  1420 

1410  GOSUB  1600:Z=Z3:IO=IO+1:GOTO  1400 

1420  Z=(Z-1)/(Z+1) :Z2=Z*Z:Z3=Z 

1430  13=9: GOSUB  1640:Z=Z+Z 

1440  IF  10=0  THEM  1460 

1450  FOR  P=l  TO  IO:Z=Z+Z:NEXT 

1460  Z2=Z*I2:RETURN 

1470  ' 

1480  '***EXPONENTIAL*** 

1490  ' 

1500  Z2=EXP(ABS(Z)) :I0=0 

1510  IF  Z*Z<.004  THEN  1530 

1520  Z=Z/2:I0=I0+l:GOTO  1510 

1530  Z2=l:FOR  12=8  TO  1  STEP  -1 : Z2=Z2*Z/I2+1 : NEXT 

1540  IF  10=0  THEN  1560 
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1550  FOR  P=l  TO  I0:Z2=Z2*Z2:NEXT 

1560  RETURN 

1570  ' 

1580  '***SQUARE  ROOTS** 

1590  ' 

1600  Z3=SQR(Z) :Z3=(Z3+Z/Z3)/2:Z3=(Z3+Z/Z3)/2:RETURN 

1610  ' 

1620  '***POWER  EXPANSION  SUBROUTINE*** 

1630  ' 

1640  Z2=Z*Z:FOR  P=3  TO  ABS(I3)  STEP  2 : Z3=SGN( 13 )*Z3*Z2 

1650  Z=Z+Z3/P: NEXT: RETURN 

1660  END 


EXPEQNS/BAS 

10  'PROGRAM  BY  R.  PERNANDEZ-PRINI  &  S.  DAVID  KLEIN,  DEP 

T.  OF  CHEM. ,  UNIV.  OF  FLA.,  GAINESVILLE,  FL   32601   USA 

20  CLS 

30  PRINT"THIS  IS  THE  EXPANDED  PITTS  EQUATION  PROGRAM." 

40  PRINT :PRINT"'RESULTS  OUTPUT  TO  PRINTER." 

50  DIM  EC(25),SL(25),G(25),CL(25),CI(25),AS(25),DL(25), 

DF(25),DS(25),AX(25) 

60  INPUT"INPUT  DC,T(DEG-K),V(POISE)";DC,T,V 

70  INPUT" INPUT  N,A?PROX  LAM-0 , APPRO X  KA";N,0L,X1 

80  PRINT"NOW  INPUT  THE  CONC.-COND.  PAIRS." 

90  FOR  1=1  TO  N 

100  INPUT  EC(I) ,EL(I) 

110  EL(I)=0.999505-EL(I) 

120  NEXT  I 

130  LPRINT"EXPANDED  EQUATION" 

140  INPUT"INPUT  AO  IN  ANGSTROMS" ;A0 

150  INPUT"INPUT  AO  INCREMENT,  # • OF  ITERATIONS" ; AI , Q 

160  INPUT"FUOSS-HSIA  OR  PITTS(F  OR  P) ,  WHICH  TREATMENT" 

•  A$ 

170  LPRINT"N="  ;N;If;   DC  =  ";DC 

180  LPRINT"T(K)=" ;T;";   V=";V  :  LPRIMT 

190  FOR  1=1  TO  N 

200  LPRINTUSING"EC  =  //#.  //#*#[[[[    EL  =  ##// .  ####"  ;  EC  ( I )  ,  EL( 

I) 

210  NEXT  I 

220  M=M+1 

230  LPRINT:LPRINT"A0=" ;A0 

240  IF  A$="P"  THEN  LPRINT"PITTS  TREATMENT" : GOTO  260 

250  IF  A$="F"  THEN  LPRIHT"FUOSS-HSIA  TREATMENT"  ELSE  CL 

S:PRINTCHR$( 23) "WHICH  TREATMENT  DO  YOU  WANT,  F(UOSS-HSI 

A)  OR  P(ITTS)":FOR  1=1  TO  1000 : NEXT: CLS : GOTO  160 

260  DT=DC*T 

270  A=82.0460E4/(SQR(DT[3)) 

280  B=82.487/(V*SCR(DT) ) 
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290  A1=1.8246e6/(SQR(DT[3)) 

300  BB=5.029160E1*A0/SQR(DT) 

310  S=A5OL+B 

320  FOR  1=1  TO  N 

330  AS(I)=EL(I)/(OL-S*SQR(EC(I)*EL(I)/OL)) 

340  NEXT  I 

350  E1=2.94257S12/(DT[3) 

360  E2  =  4.3324iJE7/(V*(DT[2)) 

370  E=E1*0L-E2 

380  BA=l6.7099E4/DT 

390  K=50.2916/SQR(DT) 

400  BZ=BA/A0 

410  IF  A$  =  "F"  THEN  GOTO  430 

420  Jl=2*El*OL*(LOG(K*A0)+2/BZ+1.77l8)+B*K*A0+2*E2*(0.0 

1387-L0G(K*A0))  :  GOTO  440  'THIS  IS  PITTS  Jl 

430  Jl=El*OL*(1.8l47+2*LOG(K*A0)+2s(2/BZ+2/(BZ[2)-l/(BZ 

[3)))-E2*(1.5337+4/(3sBZ)+2sLOG(K*A0))+A*B+B*K*A0  :  GOT 

0  450  'THIS  IS  FU0S3-HSIA  Jl 

440  J2=-(4*K*BA*El*OL*(1.2929/(BZ[2)+1.5732/BZ)+B*( (K*A 

0)[2)+1.40725-(2*E2*K*A0))  :  GOTO  460  'THIS  IS  PITTS  J2 

450  J2=-(El^OL*K*A0s(0.6094+4.4748/BZ+3.8284/(BZ[2))-E2 

*K*A0*(1.5405+2.276l/BZ)+El*B»(2*(2/BZ+2/(BZ[2)-l/(BZ[3 

) )-1.9384)-B*E2*(4*l/(3*BZ)-2.2x94)/OL+A*B«K*A0+B*(K*A0 

)[2)  'THIS  IS  FUOSS-HSIA  J2 

460  S=A*OL+B 

470  FOR  1=1  TO  N 

480  CI(I)=AS(I)*EC(I) 

490  AS(I)=EL(I)/(OL-S*SQR(CI(I))+E*CI(I)*LOG(CI(I) )+Jl* 

CI(I)+J2*SQR(CI(I)[3)) 
500  CI(I)=AS(I)*EC(I) 

510  G(I)=10[(-A1*SCR(CI(I))/(1+EB*SQR(CI(I)))) 

520  DL(I)=EL(I)-OL+S*SQR(CI(I))-E*CI(I)*LOG(CI(I))-Jl*C 

I(I)-J2*(SQR((CI(I))[3))+K1*G(I)[2*EL(I)*CI(I) 

530  NEXT 

540  S1=0:S2=0:S3=0:S4=0 

550  FOR  1=1  TO  N 

560  S1=S1-EL(I)*CI(I)*C(I)[2 

570  S2=S2+EL(I)[2*CI(I)[2*G(I)[4 

580  S3=33+DL(I) 

590  S4=S4-DL(I)*EL(I)*G(I)[2*CI(I) 

600  NEXT 

610  DZ=N*S2-S1[2 

620  LI=(S3*S2-S1*S4)/DZ 

630  KI=(N*S4-S1*S3)/DZ 

640  OL=OL+LI 

650  K1=K1+KI 

660  PRINTUSING"OL=###.###    LI  =  //#//  .#//#";  OL  , LI 

670  PRINTUSING"KA=##.  ####[[[[    KI  =  ##  .####[[[["  ;K1 ,  KI 

680  IF  ABS(LI)>0.001  THEN  GOTO  710 

690  IF  ABS(KI/K1)> 0.001  THEN  GOTO  710 

700  GOTO  720 

710  GOTO  370 

720  FOR  1=1  TO  N 

730  AX( I) = ( -1+SQR( 1+4*K1*G( I) *G( I) *EC ( I) ) )/ ( 2*K1*G( I ) *G 

(I)*EC(I)):  IF   :  '  AX(I)-AS(I))>0.0005  THEN  AS(I)=AX(I) 

:CI(I)=EC(I)*AS(I) :G(I)=10[(-A1*SQR(CI(I))/(1+BB*SQR(CI 

(I)))) :GOTC730 

740  CL(I)=A3(I)*(0L-S*SQR(CI(I))+E*CI(I)*L0G(CI(I))+J1* 

CI(I)+J2*SQR((CI(I))[3)) 
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730  AX(I)=(-1+SQR(1+4*K1*G(I)*G(I)*EC(I)))/(2*K1*G(I)*G 

(I)*EC(I)):  IF  A3S(AX(I)-AS(I) )>0.0005  THEM  AS(I)=AX(I) 

:CI(I)=EC(I)*AS(I) :G(I)=10[(-A1*5QR(CI(I))/(1+BB*3QR(CI 

(I) ) ) ) :GOTO730 

740  CL(I)=AS(I)*(OL-S*SQR(CI(I) ) +E*CI(  I) *LOG(CI(I) ) +J1* 

CI(I)+J2*SQR((CI(I))[3)) 

750  DF(I)=EL(I)-CL(I) 

760  DS(I)=DF(I)[2 

770  NEXT  I 

780  SD=0 

790  FOR  1=1  TO  N 

800  SD=SD+DS(I) 

810  NEXT  I 

820  SE=SQR(SD/(N-2)) 

830  SL=SE*SQR(S2/DZ)  :  SK=SE*SQR(N/DZ ) 

840  LPRINTUSING"LAM-0=###.###   SD  LAM-0=## .###"; OL , SL 

850  LPRINTUSING"KA=##. ####[[[[  SD  KA=##. ###[[[[" ;K1 ,SK 

860  LPRINTUSING"SD  LAM=# ##.'###"  ;SE 

870  IF  M<Q  THEN  A0=A0+AI  :  GOTO  220 

880  INPUT"DO  YOU  WISH  TO  CHANGE  TREATMENTS (Y/N) "; B$ : IF 

B$="Y"  THEN  GOTO  900 

890  GOTO  920 

900  IF  A$="P"THEN  A$="F" : A0=A0-(Q-1) *AI :M=0 : GOTO220 

910  IF  A$="F"THEN  A$=" P" : A0=A0-( Q-l ) * AI : M=0 : GOTO220 

920  END 


SHEDKAY/BAS 

10  'PROGRAM.  BY  S.  DAVID  KLEIN,  DEPT.  OF  CHEM.  ,  UNIV.  OF 

FLA.,  GAINESVILLE,  FL  32601  USA 
20  CLS 

30  PRINT"***  THIS  IS  THE  SHEDLOVSKY-KAY  EXTRAPOLATION  P 
ROGRAM  ***" : PRINT 
40  DEFDBLA-0,R-Z 

50  PRINT"ENTER  ALL  NUMBERS  IN  DOUBLE  PRECISION !": PRINT 
60  INPUT  "INPUT  CELL  CONSTANT , SOLVENT  RESIS";CC,SR 
70  INPUT"INPUT  DC,T,V" ;DC,T,V 
80  Z=DC*T:GOSUB  990:DB=Z3 
90  AL=8.203D5/( (DC*T)*DB) 
100  BE=82.43D0/(V*DB) 
110  IIIPUT"INPUT  N";N 
120  F0RQ=1T0N 

130  INPUT"IMPUT  C,RES";EC(0) ,R(Q) 
140  L(C)=1000*CC/R(0) 

150  EL(Q)=100O*CC*((l/R(Q) ) -( 1/SR) )/EC (0) 
160  NEXT 

170  L0=1000D0*CC/SR 
180  A=3.649D6/((DC*T)*DB) 
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190  INPUT"INPUT  LAH0";OL 

200  Z=OL:GOSU3990:OP=Z3 

210  F0RQ=1T0N 

220  Z=EC(Q)*EL(Q) :G0SU3  990:EK(Q)=Z3 

230  Z(Q)=(AL*OL+BE)/(OL*OP)*EK(Q) 

240  Z=1+(Z(Q)/2D0)*(Z(Q)/2D0) : GOSUB990 : SP( 0) =Z3 

250  3(Q)=(Z(Q)/2D0+SP(Q) ) * ( Z (0) /2D0+SP( Q) ) 

260  X(Q)=EL(Q)*S(Q)/OL 

270  NEXT 

280    PORQ=lTOri 

290  Z=10D0:GOSUB  820:LT=Z2 

300  Z=EC(Q)*X(Q) :GOSUB990:PM(Q)=Z3"(-A) 

310  Z=FM(Q)*LT:GOSUB910:P(Q)=Z2 

320  Z=F(Q) :GOSUB990:P(Q)=Z3 

330  K(Q)=EC(Q)*X(Q)*X(Q)*F(Q)*F(Q)/(1D0-X(Q)) 

340  Z=K(Q) :GOSUB990:KN(Q)=Z3 

350  Z=EC(Q) :GOSUB990:SE(Q)=Z3 

360  Z=(1-((L(Q)/EC(Q))*S(Q)/0L)) : GOSUB990 :0G( Q) =Z3 

370  L(Q)=OL*KN(Q)*SE(Q)/(S(Q)sF(Q))*0G(Q)+L0 

380  NEXT 

390  FORO=lTON 

400  IA(Q)=SE(Q)/(S(Q)*F(Q))*OG(Q) 

410  IB(Q)=L(Q) 

4  20  NEXT 

430  GOSUB  590 

440  PRINT"SD  OF  1ST  APPROX  LST  SQ=" ;SD 

450  L0=IH" 

460  F0RQ=1T0N 

470  Z=(l-((L(Q)/EC(Q))*S(Q)/OL)*(1.0D0-L0/L(Q))) :GOSUB9 

90:LK(Q)=Z3 

480  L(Q)=OL*KN(Q)*SE(Q)/(S(Q)*F(Q) )*LX(Q)+L0 

490  NEXT 

500  FGRQ=1T0N 

510  IA(Q)=SE(Q)/(S(Q)*F(Q))*LK(Q) 

520  IB(Q)=L(0) 

530  NEXT 

540  GOSUB  590 

550  PRINT"SD  OF  2ND  LST  SQR=";SD 

560  K=(SL/OL)*(SL/OL) 

570  PRINT "K=";K 

580  END 

590  'LEAST  SQUARE  SUBROUTINE 

600  'X=IA    Y=IB 

610  IE=0:IK=0:IG=0:IH=0:IJ=0 

620  F0RQ=1T0N 

630  IC(Q)=IA(Q)*IA(Q) 

640  ID(Q)=IB(Q)*IB(Q) 

650  IE=IE+IA(Q) 

660  IK=IK+IB(Q) 

670  IG=IG+IC(0) 

680  IH=IH+ID(Q) 

690  II(Q)=IA(Q)*IB(Q) 

700  IJ=IJ+II(Q) 

710  NEXT 

720  UZ=IE/N:VZ=IK/N 

730  SX=( :  :  N)-(UZ*UZ) 

740  SY=(IH/N)-(VZ*VZ) 

750  SL=((IJ/N)-(UZ*VZ))/SX 

760  IN=VZ-(SL*UZ) 
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730  SX=(IG/N)-(UZ*UZ) 

740  SY=(IH/N)-(VZ*VZ) 

750  SL=((IJ/N)-(UZ*VZ))/SX 

760  IN=VZ-(SL*UZ) 

770  VA=((IH-(VZ*IK))-((SL*SL)*(IG-(UZ*IE))))/(N-2) 

780  SD=VA/(IG-(UZ*IE)) 

790  RETURN 

800  END 

810  ****  NATURAL  LOGARITHM  *** 

820  Z2=LOG(Z) :I0=0:I2!=SGN(Z2) :IF  I2!<0  Z=1D0/Z 

830  IF  Z<1.065  THEN  850 

840  GOSUB  990:Z=Z3:I0=I0+lD0:GOTO830 

850  Z=(Z-1D0)/(Z+1D0) :Z2=Z*Z:Z3=Z 

860  I3=9:GOSUE  1010:Z=Z+Z 

870  IF  10=0   THEN  890 

880  FORP=lTOI0:Z=Z+Z:NEXT 

890  Z2=Z*I2! : RETURN 

900  '***  EXPONENTIAL  *** 

910  Z2=EXP(ABS(Z)) :I0=0 

920  IF  Z*Z<.004  THEN  940 

930  Z=Z/2D0:I0=I0+lD0:GOTO920 

940  Z2=l:FORI2!=8TOlSTEP-l:Z2=Z2*Z/I2!+lD0:NEXT 

950  IF  10=0  THEN  970 

960  FORP=1TOIO:Z2=Z2*Z2:NEXT 

970  RETURN 

980  '***  SQUARE  ROOT  »** 

990  Z3=SQR(Z) :Z3=(Z3+Z/Z3)/2D0:Z3=(Z3+Z/Z3)/2D0:RETURN 

1000  '  ***  POWER  EXPANSION  SUBROUTINE  *** 

1010  Z2=Z*Z:FORP=3TOABS(I3)STEP2:Z3=SGN(I3)*Z3*Z2 

1020  Z=Z+Z3/P:NEXT:RETURH 

1030  END 


FREQEXTR/BAS 

10  DEFDBL  A-H,J-Z 

20  CLS 

30  '  PROGRAM  BY  S.  DAVID  KLEIN,  UNIV.  OF  FLA.,  DEPT.  OF 

CHEM.,  GAINESVILLE,  FL   32611  USA 
40  ZZ=0 

50  PRINT"THIS  IS  THE  FREQUENCY  EXTRAPOLATION  AND  EQUIV. 
COND.  PROGRAM" 
60  PRINT"INPUT  FREQUENCY, RESISTANCE  PAIRS  (FREQ.  IN  KHZ 

)" 

70  ZZ=ZZ+1 

80  FOR  1=1  TO  3 

90  INPUT  F(I) ,R(I) 

100  FR( I )=?( I )*2*3- 141592654*1000 

110  F2(I)=FR(I)[2 

120  NEXT  I 

130  Z  =  5.00003*(R(2)-R(1))/(R(3)-R(D) 
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FREQEXTR/BAS 

10  DEFD3L  A-H,J-Z 

20  CLS 

30  '  PROGRAM  EY  S.  DAVID  KLEIN,  UNIV.  OF  FLA.,  DEPT.  OF 

CHEM.,  GAINESVILLE,  FL   32611  USA 
40  ZZ=0 

50  PRINT"THIS  IS  THE  FREQUENCY  EXTRAPOLATION  AND  EQUIV. 
COND.  PROGRAM" 
60  PRINT" INPUT  FREQUENCY, RESISTANCE  PAIRS  (FREQ.  IN  KHZ 

)" 

70  ZZ=ZZ+1 

80  FOR  1=1  TO  3 

90  INPUT  F(I) ,R(I) 

100  PR(I)=P( I) *2*3. 141592654*1000 

110  F2(I)=PR(I)[2 

120  NEXT  I 

130  Z  =  5.  00003s  (R(  2)  -R(1))/(R(3)-R(D) 

140  A=(Z-1)/(P2(3)-(Z*P2(2))) 

150  PRINT  Z,A 

160  FOR  1=1  TO  3 

170  B(I)=1/(1+(A*P2(I)) ) 

180  NEXT  I 

190  M1=(R(1)-R(2))/(B(1)-B(2)) 

200  M2=(R(2)-R(3))/(B(2)-B(3)) 

210  M3=(Ml+M2)/2 

220  IP  M3-M1> 0.001  OR  M3-M2> 0.001  THEN  48 0 

230  PRINT  "STRAIGHT  LINE  OK" 

240  IN#=R(2)-(M3*B(2)) 

250  PRINTUSING"INTERCEPT=######.##  =R(S)";IN# 

260  PRINT"NOW  INPUT  CONCENTRATION (MOL/GM  SOL'N)TO  GET  M 

OLARITY" 

270  INPUT  C 

280  IF  ZZ>1  THEN  GOTO  310 

290  PRINT"IMPUT  THE  SOLUTION  DENSITY  SLOPE, THE  SOLVENT 

DENSITY." 

300  INPUT  SL,SD 

310  MO=CS(SL*C+SD)*1000 

320  PRINTUSING"MOLARITY=  ##.####[[[[  MOLES  PER  LITRE."; 

MO 

330  IF  ZZ>1  THEN  GOTO  360 

340  PRINT"INPUT  THE  SOLVENT  RESISTANCE." 

350  INPUT  SR 

360  RS=SR*IN#/(SR-IN#) 

370  IF  ZZ>1  THEN  GOTO  410 

380  PRINT"INPUT  THE  CELL  CONSTANT." 

390  INPUT  CC 

400  IF  B$="Y"  THEN  GOTO  60 

410  LA=1000*CC/(RS*MO) 

420  PRINTUSING"EQUIV.  COND  .=#// # // .  # # (fff    S-CM[  2/MOL"  ;  LA 
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430  PRINT 

440  INPUT"CONTIHUE  WITH  SAME  RUN  AMD  SAME  CELL";A$ 

450  IP  A$="N"  THEN  PRINT"SAME  RUN,  DIFFERENT  CELL?"  ELS 

E  GOTC60 

460  INPUT  B$:IF  B$="Y"  THEN  GOTO  330 

470  END 

480  PRINT  Ml, M2, "PROBLEM  HERE!" 

490  END 
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